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Abstract of the Dissertation 

QCD resummation and heavy quark cross 

sections 

by 

Nikolaos Kidonakis 
Doctor of Philosophy 
in 
Physics 

State University of New York at Stony Brook 
1996 



In this dissertation a detailed analysis of heavy quark produc- 
tion is given with an emphasis on the resummation of soft gluon 
corrections. 

First we calculate the production cross sections for top quark 
production at the Fermilab Tevatron and for bottom quark produc- 
tion at fixed-target pp experiments, and in particular HERA-B. We 
consider both the order cross sections and the resummation of 
soft gluon corrections in all orders of QCD perturbation theory. 

Then we calculate the inclusive transverse momentum and ra- 
pidity distributions for top quark production at the Fermilab 

iii 



Tevatron and bottom quark production at HERA-B. We give both 
al and resummed results. 

Finally, we discuss the resummation of distributions that are 
singular at the elastic limit of partonic phase space (partonic thresh- 
old) in QCD hard-scattering cross sections, such as heavy quark 
production. We show how nonleading soft logarithms exponentiate 
in a manner that depends on the color structure within the un- 
derlying hard scattering. This result generalizes the resummation 
of threshold singularities for the Drell-Yan process, in which the 
hard scattering proceeds through color-singlet annihilation. We 
illustrate our results for the case of heavy quark production by 
light quark annihilation and gluon fusion, and also for light quark 
production through gluon fusion. 
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Chapter 1 



Introduction 



In this dissertation we will be concerned with the calculation of heavy 
quark cross sections and differential distributions in transverse momentum and 
rapidity. In particular we will examine top quark production at the Fermilab 
Tevatron and bottom quark production at fixed target pp experiments and 
HERA-B. 

The Standard Model of Particle Physics is at present the most successful 
model for the description of the interactions of elementary particles. In this 
model all the fundamental interactions derive from the principle of local gauge 
invariance. The Standard Model is described by the gauge group SU(3)c <8> 
SU(2)l ® U(1)y, which incorporates the Glashow-Weinberg-Salam theory of 
electroweak processes and Quantum Chromodynamics (QCD). QCD is the 
local SU(3)c gauge field theory of strong interactions and in perturbative QCD 
we make physical predictions by expanding in powers of the strong coupling 
constant a s . 

In the Standard Model there are three families of quarks and leptons. The 
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heaviest family of quarks consists of the top and the bottom quarks. The latter 
was discovered in the seventies but the top was elusive until recently due to 
its very high mass. The search for the top has gone on for more than a decade 
with experimental groups continually establishing ever higher limits for the top 
quark mass. The recent discovery of the top quark by the two experimental 
groups CDF [|TJ and DO [0 at the Fermilab Tevatron has been the biggest 
discovery in particle physics in the last decade and has given us one more 
confirmation of the Standard Model. The existence of the top was expected 
on theoretical grounds for the cancellation of anomalies. CDF announced a 
top quark mass of 176 ± 8 ± 10 GeV/c 2 and a it production cross section of 
6.8i|]4 pb. DO measured the top quark mass to be 199t$ ± 22 GeV/c 2 and 
its production cross section to be 6.4 ± 2.2 pb. 

Perturbative QCD (pQCD) is applicable only at high energies since it is 
only in this domain that a s is small according to the principle of asymptotic 
freedom ||. The expansion of a s in the renormalization scale /i is given by 

6(153 - 1971/) ln(ln(/i 2 /A 2 ))" 
(33 - 2n f ) 2 ln(/i 2 /A 2 ) J + " ' 

(1.1) 

where rif is the number of quarks with mass less than \x and A is the QCD 
scale parameter. 

In fig. 1 we plot a s as a function of scale. We see that at the scale 
relevant to top quark production (i.e. the top quark mass, which we take to 
be 175 GeV/c 2 ) the value of the coupling is about 0.1. This number is small 
enough to allow perturbative QCD calculations to be reliable. For bottom 



a s (ii 2 ,n f ) 
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(33 - 2n f ) In(/J 2 /A 2 
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Figure 1.1: The strong coupling constant a s versus scale 

quark production the relevant scale is the bottom quark mass (which we take 
to be 4.75 GeV/c 2 ). There the value of the coupling is about 0.2, which is still 
small but not as much as for the top quark case. 

The calculation of production cross sections for heavy particles in QCD 
is made by invoking the factorization theorem || and expanding the contribu- 
tions to the amplitude in powers of the coupling constant a s (fi 2 ). One has to 
perform both renormalization of ultraviolet divergences and mass factorization 
of collinear divergences. In top and bottom quark production the heavy quark- 
antiquark pairs are created in parton-parton collisions. For the Born 0(a 2 s ) 
cross section the two relevant processes are quark-antiquark annihilation 



q + q — > Q + Q, 



(1.2) 



and gluon-gluon fusion 

g + g^Q + Q. (1.3) 
In 0(a^) the relevant parton-parton processes are 

q + q^Q + Q + g, (1.4) 
3 + 2^Q + <2 + 2, (1.5) 

q + g^Q + Q + q, (1.6) 

g + 2^Q + Q + g, (1.7) 

together with the virtual corrections to (1.2) and (1.3). Recent investigations 
have shown that near threshold there are large logarithms in the perturbation 
expansion which have to be resummed to make more reliable theoretical pre- 
dictions. The standard process is fixed target Drell-Yan production, which has 
been the subject of many papers over the past few years |J. The same ideas 
on resummation were applied to the calculation of the top-quark cross section 
at the Fermilab Tevatron in || and ]7|. What is relevant in these reactions is 
the existence of a class of logarithms of the type (ln(l — z)) l /(l — z), where % 
is the order of the perturbation expansion, and where one must integrate over 
the variable z up to a limit z — 1. These terms are not actually singular at 
z = 1 due to the presence of terms in 5(1 — z). However the remainder can 
be quite large. In general one writes such terms as "plus" distributions, which 
are then convoluted with regular test functions (the parton densities). 

In perturbation theory with a hard scale we can use the standard 
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expression for the order- by order cross section in QCD, namely 
a(S, m 2 ) = / 2 dxi f 2 dx 2 fi( x i^ 2 )fj( x 2, l^ 2 )^ij(s = x x x 2 S, m 2 , fx 2 ) , 

S 11S ij 

(1.8) 

where the fi(x,fi 2 ) are the parton densities at the factorization scale fi 2 and 
the Oij are the partonic cross sections. The numerical results for the hadronic 
cross sections depend on the choice of the parton densities, which involves 
the mass factorization scale fi 2 ; the choice of the running coupling constant, 
which involves the renormalization scale (also normally chosen to be fi 2 ); and 
the choice for the actual mass of the heavy quark. In lowest order or Born 
approximation the actual numbers for the cross section show a large sensitivity 
to these parameters. In chapter 2 we will show plots of the top and bottom 
quark production cross sections in leading order (LO), i.e. 0(a 2 ), and next-to- 
leading order (NLO), i.e. 0(a^). The NLO results follow from the work of the 
two groups H and ||. However even including the NLO corrections does not 
completely fix the cross section. The sensitivity to our lack of knowledge of 
even higher terms in the QCD expansion is usually demonstrated by varying 
the scale choice up and down by factors of two. In general it is impossible to 
make more precise predictions given the absence of a calculation in next-to- 
next-to-leading order (NNLO). However in specific kinematical regions we can 
do so. 

The threshold region is one of these regions. In this region one finds that 
there are large logarithms which arise from an imperfect cancellation of the 
soft-plus-virtual (S+V) terms in the perturbation expansion. These logarithms 
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are exactly of the same type mentioned above. In || the dominant logarithms 
from initial state gluon bremsstrahlung (ISGB), which are the cause of the 
large corrections near threshold, were carefully examined. Such logarithms 
have been studied previously in Drell-Yan (DY) || production at fixed target 
energies (again near threshold) where they are responsible for correspondingly 
large corrections. The analogy between DY and heavy quark production cross 
sections was exploited in || and a formula to resum the leading logarithms in 
pQCD to all orders was proposed. Since the contributions due to these loga- 
rithms are positive (when all scales \i are set equal to the heavy quark mass 
m), the effect of summing the higher order corrections increases the top and 
bottom quark production cross sections over those predicted in 0(a^). This 
sum, which will be identified as <T res , depends on a nonperturbative parameter 
//o- The reason that a new parameter has to be introduced is that the resum- 
mation is sensitive to the scale at which pQCD breaks down. As we approach 
the threshold region other, nonperturbative, physics plays a role (higher twist, 
bound states, etc.) indicated by a dramatic increase in a s and in the resummed 
cross section. This is commonly called the effect of the infrared renormalon or 
Landau pole flO]. We choose to simply cut off the resummation at a specific 
scale /jlq where A « fi « m since it is not obvious how to incorporate 
the nonperturbative effects. Note that our resummed corrections diverge for 
small fio but this is not physical since they should be joined smoothly onto 
some nonperturbative prescription and the total cross section will be finite. 
Another way to make it finite would be to avoid the infrared renormalon by a 
specific continuation around it, i.e. the principal value resummation method 
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H|. However, at the moment our total resummed corrections depend on the 



parameter fiQ for which we can only make a rough estimate. 

We will see in chapter 2 that the gluon-gluon channel is the dominant 
channel for the production of 6-quarks near threshold in a fixed-target pp ex- 
periment. This is not the case for the production of the top quark at the Fer- 
milab Tevatron, which is a proton-antiproton collider, and where the dominant 
channel is the quark-antiquark one. That was fortunate as the exponentiation 
of the soft-plus-virtual terms in is on a much more solid footing in the qq 
channel, due to all the past work which has been done on the Drell-Yan reac- 
tion [||]. We will examine all "large" corrections near threshold, including both 
Coulomb-like and large constant terms. We will do that in chapter 2 where we 
will present all the relevant formulae at the partonic level. In addition we will 
present subleading S+V terms and discuss their contribution to the total S+V 
cross section. Chapter 2 also contains the analysis of the hadron-hadron cross 
section which is relevant for top quark production at the Fermilab Tevatron 
as well as for bottom quark production for the HERA-B experiment and for 
fixed-target pp experiments in general. We give results in LO, in NLO and 



after resummation. Most of our results in chapter 2 have appeared in [0 and 

In chapter 3 we will present results for the heavy quark differential dis- 
tributions in transverse momentum pt and rapidity Y . After showing the 
relevant formulas at the partonic level and discussing resummation we will 
give results for top and bottom quark production at the Fermilab Tevatron 
and HERA-B, respectively. We will show that resummation produces an 
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enhancement of the NLO results with little change in shape. Our results 



in this chapter have appeared in and [ |T4|j . 

In chapter 4 we will present a new resummation formula that takes into 
account subleading logarithms. A strong motivation for this research is the 
inadequacy of the leading-log approximation for the gg channel as will be 
discussed in chapter 2. We will show how nonleading soft logarithms exponen- 
tiate in a manner that depends on the color structure within the underlying 
hard scattering. We will derive the anomalous dimension matrices and ex- 
hibit the exponentiation of nonleading soft logarithms for the production of 
heavy quarks through both light quark- ant iquark annihilation and gluon-gluon 



fusion. Some of our results in chapter 4 have been presented in [|15 
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Chapter 2 

Top and bottom quark production cross 
sections 

In this chapter we present general formulas for heavy quark cross sections. 
First we discuss our results at the partonic level and then we give the calcula- 
tions for the cross sections for top quark production at the Fermilab Tevatron 
and for bottom quark production at fixed-target energies and HERA-B. We 
present both the order cross sections and the resummation of soft gluon 
corrections in all orders of QCD perturbation theory. 

2.1 Results for parton-parton reactions 

The partonic processes that we examine are 

i(h)+j(k 2 )^Q( Pl ) + Q( P2 ), (2.1.1) 

where i,j = g,q,q. The square of the parton-parton center-of-mass (cm.) 
energy is s = (ki + k^) 2 . 
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We begin with an analysis of heavy quark production in the qq channel. 
The Born cross section in this channel is given by 

4W) = f al{^)K qq NC F l -P (l + *f\ (2.1.2) 

where Cp = (N 2 — 1)/(2N) is the Casimir invariant for the fundamental 
representation of SU(N), K qg - = N~ 2 is a color average factor, \i denotes the 



renormalization scale, and (3 = yl — Am? / s. Also iV = 3 for the SU(3) color 
group in QCD. The threshold behavior (s — > Am 2 ) of this expression is given 
by 



a 



(°) (c — 
qq, thres 



s,m 2 ) = ira 2 s (/j 2 )K q qNC F -[3. (2.1.3) 

Complete analytic results are not available for the NLO cross section as some 
integrals are too complicated to do by hand. However in analytic results are 
given for the soft-plus-virtual (S+V) contributions to the cross section, and 
for the approximation to the cross section near threshold. Simple formulas 
which yield reasonable approximations to the exact 0(a^) results have been 
constructed in . From these results one can derive that the Coulomb terms 
to first order in the qq channel are given by 

«£■>(., m") = (C F - St) (2.1.4) 

in the MS scheme, where Ca = N is the Casimir invariant for the adjoint 
representation of SU(N). These terms are distinguished by their typical 
behavior near threshold which, after multiplication by the Born cross section, 
yield finite cross sections at threshold in NLO. We note that Cp— Ca/2 = — 1/6 
is negative for 577(3) and that the first-order Coulomb correction is negative 
(the interaction is repulsive). 
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In the DIS scheme, in addition to the Coulomb terms, we also have a large 
constant contribution so that the first order result near threshold is 



a\ q ; (s,m ) = (Tqq'{S,m 



2(3 V 2 ) 71 V2 3 



(2.1.5) 

We have included the constant terms to see their effect at larger values of j3. 

Since the total parton-parton cross sections only depend on the variables 
s and m 2 they can be expressed in terms of scaling functions as follows 



avis, to 2 ) = ^(r{fO,m 2 ) 

k=0 

= ^ i E(4-« s (/i 2 )) fc E/f ) ^ ln ^' (2.1.6) 

where we denote by the 0(a k s +2 ) contribution to the cross section. The 
scaling functions fij ,l \rj) depend on the scaling variable r] = s/Am 2 — 1 = 



s(3 2 /Am 2 . 



In fig. 2.1 we plot fqg (t]) for k = 0, 1 for the exact and threshold 
expressions (from |l[]) in the MS scheme. We see that the threshold Born 
approximation is excellent for small rj and reasonable for the entire range of rj 
shown. We also note that the theshold first-order approximation is good only 
very near to threshold. In fig. 2.2 we plot the corresponding functions for 
the DIS scheme. Here the first-order corrections are smaller than in the MS 
scheme. Again the threshold first-order approximation is good only very close 
to threshold. 

The analysis of the contributions to the gluon-gluon channel in NLO is 
much more complicated. First of all there are three Born diagrams each with 



0.0001 0.001 0.01 0.1 1 10 

r/ = s/ 4rri s 1 

Figure 2.1: The scaling functions /^'°' ) in the MS scheme. Plotted are / < ^'°' ) 
(exact, upper solid line at large 77; threshold approximation, upper dotted line 
at large 77), /^'°' ) (exact, lower solid line at large 77; threshold approximation, 
lower dotted line at large 77). 




Figure 2.2: Same as fig. 2.1 but now for the DIS scheme. 
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a different color structure. Therefore only few terms near threshold are pro- 
portional to the Born cross section. The exact Born term in the gg channel 
is 



47rai(fi 2 )K gg NC F - lC F 



4m 



-(1 + 
In 



6m 



+ c. 



+ 1+ 



1 5 m 2 \ „ 4m 4 
3 + 3- ) l3+ -^ 



In 



4m 2 \ 

s J 
l+£ 
1-/3 

l + P' 



P 



(2.1.7) 



where K gg = (N 2 — 1) 2 is a color average factor. The threshold behavior 



[s — > 4m 2 ) of this expression is given by 



a 



(0) 

gg, thresV°> 



s, m 



) = naj (fi 2 )K gg -NC F [4C F - C A ](5. 



(2.1.8) 



Again, the complete NLO expression for the cross section in the gg channel 
is unavailable but analytic results are given for the S+V terms in |f|. These 
were used in [|| to analyze the magnitude of the cross section near threshold. 
From the approximate expressions given in one can extract the 7r 2 terms to 
first order in the gg channel. These are 



af;\s,m 2 ) = al(fx 2 )NC K K : 



71 



99 



5 + -8 2 + 16 m< 



8 24 



,4\ x 



m m , , 

3 V~ 1( V 1-1,1 



1 + 



71 



+® 3 s CQ, ED K gg — 



1 m 6 
-- — lb— r- 



/ m 8 m 4 " 



f3 1-/3 



(2.1.9) 



where C# = (iV 2 — 1)/N = 2NC F C A - 4NC 2 F , and C QED = (iV 4 - l)/iV 2 = 
—4C F + ACaC f . These are not proportional to the Born term. The threshold 
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behavior of (2.1.9) is given by 

-~NC K 



T gg,L™( S > m ) = a sW ) K 99^ ~ 



+ a 



QED 



(2.1.10) 



which is proportional to the threshold Born term. Therefore the threshold 
approximation for the tc 2 terms in NLO can be written as 

rrW+^fo m 2 ) -rr (0) ( « m 2 ) [l I 7 I^l} ( ~NC k /2 + C qed \ 
°i»,th» [s,m ) - a gg , thres {s, m)^ + ^ y ^ _ ) , 

(2.1.11) 

or, writing the color factors in terms of N, as 

7ia s (fi 2 ) N 2 + 2 



rr (0)+(7r } fs m 2 ) - rr (0) (s m 2 ) 

gg, thres J — thresh J 



1 + 



(2.1.12) 



4/? iV(iV 2 - 2) 

In fig. 2.3 we plot the scaling functions fgg'°\f]) with k = 0, 1 in the 
MS scheme for the exact and threshold expressions (from ||). We see that 
the Born and first-order threshold approximations are good only very close to 
threshold. 

In H an approximation was given for the NLO soft-plus-virtual contribu- 
tions and the analogy with the Drell-Yan process was exploited to resum them 
to all orders of perturbation theory. The S+V approximation is adequate in 
the region 0.1 < r\ < 1 (which is the kinematical region of interest as we will see 
in sections 2.2 and 2.3) for the qq channel, but not as good for the gg channel 
in the MS scheme. Therefore we reexamined the approximate formulae given 
in m for the initial state gluon bremsstrahlung (ISGB) mechanism to see if 
there are subleading terms that will improve the S+V approximation. Let us 
see the structure of these terms. We are discussing partonic reactions of the 
type i(ki) -\- j{kz) —* Q{pi) + Q(P2) + g(k 3 ), and we introduce the kinematic 




Figure 2.3: Same as fig. 2.1 but now for f^ k ^ in the MS scheme. 
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variables t 1 = {k 2 — p 2 ) 2 — m 2 , U\ = {k\ —P2) 2 — m 2 , and s 4 = s + ti +u±. The 
variable s 4 depends on the four-momentum of the extra partons emitted in the 
reaction. We write the differential cross section in order «J(/i 2 ) as follows 



Here a small parameter A has been introduced to allow us to distinguish 
between the soft (s 4 < A) and the hard (s 4 > A) regions in phase space. 
The quantities a/(/i 2 ) contain terms involving the QCD /5-functions and color 
factors. The variables t\ and u\ are then mapped onto the variables s 4 and 
cos 9, where 9 is the parton-parton cm. scattering angle: 



2 d 2 a ( ij '(s,t 1 ,ui) 



a: 



2k-l 1 

j^q L s 4 v m z/ 




(2.1.13) 




(2.1.14) 



(2.1.15) 



After explicit integration over the angle 9, the series becomes 




x[alf(s,s„m 2 )-alf(s,0,m 2 )] 




where 



<t§\ 8 , S4 , m 2 ) = ±[(8- s 4 ) 2 - Asm 2 }'/ 2 £ 



1 



dcos9 ofAs, s 4 ,cos#) 



(2.1.17) 
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The Born approximation differential cross sections can be expressed by 

(2.1.18) 



2 d 2 a^ ) (s,t 1 ,u 1 ) B 

- — d{s + ti +Ui) (TiASjt^Ui) 



with 



and 



GqJsit^ui) = vra s (/i )K q qNC F 



t\ + u{ 2m 



+ 



(2.1.19) 



2na 2 ( f i 2 )K gg NC F 



c P -a 



tiUi 



X 



ti u\ Am s 
_+ + 



(1 



m 2 s 



The first-order approximate S+V result for the qq channel in 
scheme is 

2C F 



(2.1.20) 
the MS 



2 d 2 a^(s,t l ,u l ) 



dt\du\ 



7T 



-a s (fjL 



x 



- 21n^ + ln— 9(s A - A) 
Sa \ m z fi z 



+ I In. 

+ 



2 A A m > . 

- — + In — In — (5s 4 

Ca l n/ A n Ca , A 

s 4 - A) - — - In -^d{SA_) 



2C F s A 



2C f 



(2.1.21) 



The terms in the first pair of square brackets in (2.1.21) are the leading S+V 
terms given in [|J and those in the second pair of square brackets are subleading 



terms that we want to examine. In fig. 2.4 we plot the scaling functions /< 



(1,0) 



for the exact result, the leading approximate S+V result, and the approximate 
S+V result with both leading and subleading logarithms. The leading S+V 
result is a reasonable approximation to the exact result in our region of interest 
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Figure 2.4: The scaling functions fqq ' in the MS scheme. Plotted are the 
exact result (solid line), the leading S+V approximation (dotted line), the 
leading S+V approximation plus Coulomb terms (short-dashed line), the S+V 
approximation with leading plus subleading terms (long dashed line) , and the 
S+V approximation with leading plus subleading terms plus Coulomb terms 
(dash-dotted line). 
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0.1 < r] < 1. The addition of Coulomb terms worsens the leading S+V result. 
We also see that when we include the subleading terms our approximation 
does not improve much in the region of interest. However, when we add both 
the first order Coulomb terms and subleading terms to the leading S+V result 
we get a very good agreement with the exact result. In the DIS scheme the 
analogous results are 



da, 



qq 



dt\du-\ 



(s,ti,Mi) 



qq\"i 



2Cr 



7T 



S4 y m z 



9(s 4 -A) 



(1 2 A A m 2 \ 
+ - In 2 — + In — In — S(s 4 ) 



+ 



(2.1.22) 



In fig. 2.5 we plot the corresponding scaling functions. Here the addition 
of subleading terms worsens the leading S+V approximation. The addition 
of Coulomb terms and large constants enhances the first-order approximate 
results considerably. 

The resummation of the leading S+V terms has been given in The 
result is 



dt\du\ 



a q g(s,t l ,u lj 



df(si/m 2 , m 2 /V 
dsA 



-9{s A - A) 



(2.1.23) 




Figure 2.5: Same af fig. 2.4 but now for the DIS scheme. 
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where 



/_S4_ m_\ _ r Cp fs^ 2 \ 2 ^4 



7r Vm 2 ' / rn? 



r(l + 77) 

(2.1.24) 



The straightforward expansion of the exponential plus the change of the argu- 
ment in d s via the renormalization group equations generates the correspond- 
ing leading logarithmic terms. The scheme dependent A and a s in the above 
expression are given by 

Air 

A = 2, i.M>.,(»"V) ^ hW/Ar (2.1.25) 
in the MS scheme, and 

4tt 

A=l, =«.(«.') = A Kw , VAr (2.1.26) 

in the DIS scheme, where /3o = 11/3 Ca — 2/3 n/ is the lowest order coefficient 
of the QCD (3- function. The color factors CV,- are defined by C qq — C F and 
C gg = Ca, and ^ E is the Euler constant. The quantity i] is given by 



" = lf ln 1+ "°^rV>' (2L27) 



As the NNLO cross section is not known exactly we do not how to resum the 
subleading terms. In chapter 4, however, we will derive their exponentiation 
using a different formalism. 

Now let us see the analogous results for the gg channel in the MS scheme. 
We have 

2 d 2 a g y(s,t 1 ,u 1 ) 2C A 2 

S dt ldUl = «l) — «.(*0 
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Figure 2.6: Same as fig. 2.4 but now for fj^' ^ in the MS scheme. 
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(2.1.28) 



Again, the terms in the first pair of square brackets in (2.1.28) are the leading 
S+V terms and those in the second pair of square brackets are subleading 
terms. In fig. 2.6 we plot the scaling functions f^' ^ for the exact result, the 
leading approximate S+V result, and the approximate S+V result with both 
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leading and subleading terms. We note that the leading S+V approximate 
result is significantly smaller than the exact result and that the addition of 
subleading terms improves the approximation considerably. This is important 
since, as we will see in the next section, the gg channel is dominant for the 
production of 6-quarks at HERA-B. Also the addition of Coulomb terms fur- 
ther improves the approximation. The resummation of the leading S+V terms 
for the gg channel has also been given in Q. The result is 

72 - res ^ + - ^ r df(s 4 /m 2 ,m 2 //i 2 ) 



' dt ldUl = 



-0(s 4 - A) 



dS4 

2 \ 

TO \ 

(2.1.29) 



v m 2 /j, 2 

where the function / is the same as for the qq channel in the MS scheme 
(with the substitution Cp — > Ca)- Again, we do not know how to resum the 
subleading logs, but in chapter 4 we will derive their exponentiation using a 
different formalism. 

After we map t\ and u\ onto s 4 and cos# and we integrate over 6, as we 
saw earlier, the resummed cross section for either channel becomes 

. , rs-2msV 2 y p f S 4 m 2 \ d _ (0 ) . „ , 

a ij {s,m) = - ds A } — , — — er>/(s, s 4 , m ). (2.1.30) 

Note that we now have cut off the lower limit of the s 4 integration at s 4 = s 
because a s diverges as s 4 — > 0. This parameter so must satisfy the conditions 
< So < s — 2ms 1 / 2 and Sq/itl 2 « 1. It is convenient to rewrite So in terms 
of the scale /i as 

/ 2\ 3 / 2 

■^=(^1 (MS scheme); (2.1.31) 
mr V/i 1 ) 



27 



m 2 /i 2 



(DIS scheme). 



(2.1.32) 



Here /j,q is a nonperturbative parameter |4| satisfying A 2 << /j 2 , << n 2 . 

2.2 Results for top quark production at the 
Fermilab Tevatron 

In this section we examine the production of top quarks at the Fermilab 
Tevatron. Following the notation in |3| the total hadron-hadron cross section 
in order a k s is 



where S is the square of the hadron-hadron cm. energy and i,j run over q, q 
and g. The parton flux 3>y(T, yU 2 ) is defined via 



and Hij is a product of the scale-dependent parton distribution functions 
f^x,^ 2 ), where h stands for the hadron which is the source of the parton i 



The mass factorization scale /x is chosen to be identical with the renormaliza- 
tion scale in the running coupling constant. 

In the case of the all-order resummed expression the lower boundary 
in (2.2.1) has to be modified according to the condition sq < s — 2ms 1 / 2 . 




(2.2.1) 




(2.2.2) 



H ij (xi,x 2 ,fJ 2 ) = f^ 1 {x 1 ,ii 2 )fj 2 (x 2 ,fJ l 2 ) • 



(2.2.3) 
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Resumming the soft gluon contributions to all orders we obtain 

a% s (S,m 2 ) I" dr^^^a^rS^m 2 ^ 2 ), (2.2.4) 

J To 

where is given in (2.1.30) and 

ro= [^+K + *o) 1/2 ] 2 . (2 . 2 . 5) 

Here so (or equivalently //q) is the non-perturbative parameter used to cut off 
the resummation since the resummed corrections diverge for small Sq. 

We now specialize to top quark production at the Fermilab Tevatron where 
= 1.8 TeV. Taking the top quark mass as m t = 175 GeV/c 2 then the 
ratio of m t /\/~S ~ 0.1. If we choose the renormalization scale in the running 
coupling constant as m t then a s (m 2 ) ^0.1 so a s (m 2 ) \n(\/S/m t ) ~ 0.2. This 
number is small enough that we expect a reasonably convergent perturbation 
series. In the presentation of our results for the exact, approximate, and 
resummed hadronic cross sections we use the MRSD_' parametrization for the 
parton distributions ||. Note that the hadronic results only involve partonic 
distribution functions at moderate and large x, where there is little difference 
between the various sets of parton densities. We have used the MRSD_' set 
34 as given in PDFLIB || in the DIS scheme with the number of active light 
flavors rif = 5 and the QCD scale A 5 = 0.1559 GeV. We have used the two- loop 
corrected running coupling constant as given by PDFLIB. 

First, we discuss the NLO contributions to top quark production at the 
Tevatron using the results in 0] and Jl], ||. Except when explicitly stated 
otherwise we will take the factorization scale fi = m t . In fig. 2.7 we show the 
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Figure 2.7: Fractional contributions of the qq (DIS scheme, long-dashed line) 
and gg (MS scheme, short-dashed line) channels to the total 0(a z s ) top quark 
production cross section at the Fermilab Tevatron as a function of the top 
quark mass. 



30 



2 
1.9 
1.8 
1.7 
1.6 

o 

| 1.5 
1.4 
1.3 
1.2 
1.1 

1 

150 160 170 180 190 200 210 220 

Top quark mass (GeV/c 2 ) 

Figure 2.8: The K factors as a function of top quark mass for top quark pro- 
duction at the Fermilab Tevatron for the qq channel (DIS scheme, long-dashed 
line), the gg channel (MS scheme, short-dashed line), and their sum (solid 
line). 

relative contributions to the NLO cross section of the qq channel in the DIS 
scheme and the gg channel in the MS scheme as a function of the top quark 
mass. We see that the qq contribution is the dominant one and it is about 
90% of the total NLO cross section for m t = 175 GeV/c 2 . The gg contribution 
is smaller and makes up the rest of the cross section. The contributions of the 
gq and the gq channels are negligible and are not shown. 
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In fig. 2.8 we show the K factors for the qq and gg channels and for 
their sum as a function of the top quark mass. The K factor is defined by 
K = (cr(°) + |exact)/ cr< ~°' ) ) where is the Born term and | cxact is the 
exact first order correction. We notice that the K factor is quite large for the 
gg channel, which means that higher order effects are more important for this 
channel than for qq. However, since the qq channel is dominant, the K factor 
for the sum of the two channels is only slightly larger than that for qq. These 
large corrections come predominantly from the threshold region for top quark 
production where it has been shown that initial state gluon bremsstrahlung 
(ISGB) is responsible for the large corrections at NLO [Q. This can easily be 
seen in fig. 2.9 where the Born term and the O(af) cross section are plotted 
as functions of ?7 cu t for the qq and gg channels, where r\ = (s — Am 2 ) / km 2 is the 
variable into which we have incorporated the cut in our programs for the cross 
sections. As we increase 77 cut the cross sections increase. The cross sections 
rise sharply for values of ?7 cut between 0.1 and 1 and they reach a plateau at 
higher values of r] cut indicating that the threshold region is very important and 
that the region where s » Am 2 only makes a small contribution to the cross 
sections. This is the reason why we stressed in section 2.1 that our region of 
interest for comparison of the various approximations at the partonic level was 
0.1 < 7] < 1. Note that in the last figure as well as throughout the rest of this 
section we are assuming that the top quark mass is mt = 175 GeV/c 2 . 

Next, we discuss the scale dependence of our NLO results. In fig. 2.10 
we show the 0(a^) cross section as a function of the factorization scale for the 
qq and gg channels. As the scale decreases, the Born cross section increases 



32 




Figure 2.9: Cross sections for top quark production at the Fermilab Tevatron 
versus r] cut with m t = 175 GeV/c 2 for the qq channel in the DIS scheme and 
the gg channel in the MS scheme. Plotted are the Born term (qq, upper solid 
line; gg, lower solid line) and the O(af) cross section (qq, upper dashed line; 
gg, lower dashed line). 
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Figure 2.10: The scale dependence of the cross section for top quark production 
at the Fermilab Tevatron with m t = 175 GeV/c 2 . Plotted are the 0(af) cross 
section for the qq channel (upper dashed line) and for the gg channel (lower 
dashed line). 
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without bound but the exact first order correction decreases faster so that the 
NLO cross section peaks at a scale close to half the mass of the top quark and 
then decreases for smaller values of the scale. For both the qq and gg channels 
the NLO cross section is relatively flat. Thus the variation in the NLO cross 
section for scales between m f /2 and 2m t is small. In fig. 2.11 we examine 
the Ho dependence of the resummed cross section for the qq and gg channels. 
We also show, for comparison, the fi dependence of + | app | app 
where we have imposed the same cut on the phase space of S4 (54 > sq) 
as for the resummed cross section. Here | app and <r^ | app denote the 
approximate first and second order corrections, respectively, where only soft 
gluon contributions are taken into account. The effect of the resummation 
shows in the difference between the two curves for each channel. At small 
/io, cr res diverges, signalling the presence of the infrared renormalon. There 
is a region for each channel where the higher-order terms are numerically 
important. At high values of /io the two lines for each channel are practically 
the same. For the qq channel in the DIS scheme the resummation is successful 
in the sense that there is a relatively large region of /io where resummation is 
well behaved before we encounter the divergence. For the gg channel, however, 
the situation is not as good. From these curves we choose what we think 
are reasonable values for /z - We choose /i = 8.75 GeV/c 2 (0.05 m t ) and 
17.5 GeV/c 2 (0.1 m t ) for the qq channel and /io = 35 GeV/c 2 (0.2 m t ) and 
43.75 GeV/c 2 (0.25 m t ) for the gg channel, which are the choices made in || 
corresponding to upper and lower values for the cross section, respectively. 
Note that /zq need not be the same in the qq and gg reactions because the 
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Figure 2.11: The /x dependence of the resummed cross section for top quark 
production at the Fermilab Tevatron with m t = 175 GeV/c 2 for the qq channel 
in the DIS scheme and the gg channel in the MS scheme. Plotted are <r^| s 
(upper solid line at high /i ), <r^ s (lower solid line at high /x ), and the sums 
a m + a S lap P +°"S } Upp (upper dotted line) and afj + <jW | app +a$ | app 
(lower dotted line). 
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convergence properties of the QCD perturbation series could be different in 
these channels and moreover depend on the factorization scheme. 

Since we know the exact 0(a^) result, we can make an even better esti- 
mate by calculating the perturbation theory improved cross sections defined 
by 

^-imp _ ^-res , (1) I (1) I (9 9 f{\ 

° H ~ H ~r ° H lexact ° H |app 5 \Z.4.V) 

to exploit the fact that | cxact is known and | app is included in a T £ s . 

The value of the NLO cross section for the production of a top quark 
with a mass of 175 GeV/c 2 at the Fermilab Tevatron with = 1.8 TeV is 
4.8 pb. The upper and lower values of the resummed cross section are 5.6 pb 
and 4.9 pb, respectively. The upper and lower values of the improved cross 
section are 5.8 pb and 5.1 pb, respectively. Finally, in fig. 2.12 we show the 
dependence of the top quark production cross section at the Fermilab Tevatron 
on the top quark mass. Several theoretical curves [8-10] are compared with 
recent experimental results from the DO Collaboration. 

2.3 Results for bottom quark production at 
fixed-target pp experiments and HERA- 
B 

In this section we examine the production of 6-quarks in a situation where 
the presence of large logarithms is of importance, namely in a fixed-target ex- 
periment to be performed in the HERA ring at DESY. This actual experiment 
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Figure 2.12: The dependence on the top mass of the cross section for top quark 
production at the Fermilab Tevatron. 
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has the name HERA-B (TT|, [12] and involves colliding the circulating proton 
beam against a stationary copper wire in the beam pipe. The nominal beam 
energy of the protons is 820 GeV, so that the square root of the cm. energy 
is v^S* = 39.2 GeV. Taking the 6-quark mass as rrib = 4.75 GeV/c 2 then the 
ratio of rrib/y/S ~ 1/8. If we choose the renormalization scale in the running 
coupling constant as m& then a s (ml) ~ 0.2 so a s (ml) \n(\/S / m b ) « 0.4. This 
number is small enough that we expect a reasonably convergent perturbation 



series. 



In the presentation of our results for the exact, approximate, and re- 
summed hadronic cross sections we use again the same set of MRSD_' parton 
distributions as for top quark production. In this case the number of active 
light flavors is rif = 4. 

First, we discuss the NLO contributions to bottom quark production at 
HERA-B. Except when explicitly stated otherwise we will take the factoriza- 
tion scale fi = mb where m& is the 6-quark mass. In fig. 2.13 we show the 
relative contributions of the qq channel in the DIS scheme and the gg channel 
in the MS scheme as a function of the bottom quark mass. We see that the 
gg contribution is the dominant one, lying between 70% and 80% of the total 
NLO cross section for the range of bottom mass values given. The qq contribu- 
tion is smaller and makes up most of the remaining cross section. The relative 
contributions of the gq and the gq channels in the DIS scheme are negative 
and very small and they are also shown in the plot. The situation here is the 
reverse of what we saw in the previous section for top quark production at the 
Fermilab Tevatron where qq is the dominant channel with gg making up the 
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Figure 2.13: Fractional contributions of the gg (MS scheme, short-dashed 
line), qq (DIS scheme, long-dashed line), qg (DIS scheme, lower dotted line), 
and qg (DIS scheme, upper dotted line) channels to the total 0(a 3 s ) 6-quark 
production cross section at HERA-B as a function of 6-quark mass . 
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remainder of the cross section, and gq and gq making an even smaller relative 
contribution than is the case for bottom quark production at HERA-B. The 
reason for this difference between top quark and bottom quark production is 
that the Tevatron is a pp collider while HERA-B is a fixed-target pp exper- 
iment. Thus, the parton densities involved are different and since sea quark 
densities are much smaller than valence quark densities, the qq contribution to 
the hadronic cross section diminishes for a fixed-target pp experiment relative 
to a pp collider for the same partonic cross section. 

In fig. 2.14 we show the K factors for the qq and gg channels and for their 
sum as a function of the bottom quark mass. We notice that the K factor is 
quite large for the gg channel, which means that higher order effects are more 
important for this channel than for qq. Since gg is the more important channel 
numerically, the K factor for the sum of the two channels is also quite large. 
We also show the K factor for the total which is slightly lower since we are 
also taking into account the negative contributions of the qg and qg channels. 

As in the case of top quark production these large corrections come pre- 
dominantly from the threshold region. This can easily be seen in fig. 2.15 
where the Born term and the O(a^) cross section are plotted as functions of 
r] cut for the qq and gg channels. As we increase ?7 cut the cross sections increase. 
As for the top, the cross sections rise sharply for values of r] cut between 0.1 
and 1 and they reach a plateau at higher values of r] cut indicating that the 
threshold region is very important and that the region where s » 4m 2 only 
makes a small contribution to the cross sections. Note that in the last figure 
as well as throughout the rest of this section we are assuming that the bottom 
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Figure 2.14: The K factors as a function of 6-quark mass for 6-quark production 
at HERA-B for the gg channel (MS scheme, short-dashed line), the qq channel 
(DIS scheme, long-dashed line), the sum of the gg and qq channels (dotted 
line), and the sum of all channels (solid line). 
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Figure 2.15: Cross sections for 6-quark production at HERA-B versus i] cut with 
rrib = 4.75 GeV/c 2 for the qq channel in the DIS scheme and the gg channel in 
the MS scheme. Plotted are the Born term (gg, upper solid line at high r] cut ; 
qq, lower solid line at high ^ cut ) and the 0(af) cross section (gg, upper dashed 
line; qq, lower dashed line). 
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Figure 2.16: The scale dependence of the cross section for 6-quark production 
at HERA-B with m b = 4.75 GeV/c 2 . Plotted are the 0(a^) cross section for 
the qq channel (lower dashed line) and for the gg channel (upper dashed line). 

quark mass is = 4.75 GeV/c 2 . 



Next, we discuss the scale dependence of our NLO results. In fig. 2.16 we 
show the 0(al) cross section as a function of the factorization scale for the qq 
and gg channels. We see that the NLO cross section peaks at a scale close to 
half the mass of the bottom quark and then decreases for smaller values of the 
scale. For the qq channel the NLO cross section is relatively flat. The situation 
is much worse for the gg channel, however, since the peak is very sharp and 
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Figure 2.17: The total Born (dotted line) and 0(a 3 s ) (// = m b solid line, 
\i = m b /2 upper dashed line, and /i = 2m b lower dashed line) 6-quark produc- 
tion cross sections at fixed-target pp experiments versus beam momentum for 
m b = 4.75 GeV/c 2 . 

the scale dependence is much greater. Since the gg channel dominates, this 
large scale dependence is also reflected in the total cross section. Thus the 
variation in the NLO cross section for scales between m b /2 and 2m b is large. 
For comparison we note that, as we saw in the previous section, the scale 
dependence for top quark production at the Fermilab Tevatron for m t = 175 
GeV /c 2 is much smaller. 
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In fig. 2.17 we plot the Born contribution for [i = m b and the NLO cross 
section for ji = rrib/2, nib, and 2nib, as a function of the beam momentum for 
6-quark production at fixed-target pp experiments. The big width of the band 
reflects the large scale dependence that we discussed above. We see that the 
NLO cross section is almost twice as big as the Born term for the whole range 
of beam momenta that we are showing, and in particular for 820 GeV/c which 
is the value of the beam momentum at HERA-B. We also give the NLO results 
for the individual channels in fig. 2.18 for [i = nib. 

In fig. 2.19 we examine the /i dependence of the resummed cross section. 
We also show, for comparison, the /i dependence of + crW | app | app 
with the same cut s 4 > s . The effect of the resummation shows in the 
difference between the two curves for each channel. At small /i , o" rcs diverges 
signalling the divergence of the running coupling constant. There is a region 
for each channel where the higher-order terms are numerically important. At 
large values of /x the two lines for each channel are practically the same. For 
the qq channel in the DIS scheme the resummation is successful in the sense 
that there is a relatively large region of /io where resummation is well behaved 
before we encounter the divergence. This region is reduced for the qq channel 
in the MS scheme. For the gg channel, however, this region is even smaller. 

From these curves we choose what we think are reasonable values for /i . 
We choose /i = 0.6 GeV/c 2 for the qq channel in the DIS scheme (^o/m b « 
13%) and /i = 1.7 GeV/c 2 for the gg channel (/i /nib ~ 36%). The values we 
chose for the qq and gg channels are such that the resummed cross sections are 
slightly larger than the sums a^+a^ | app | app . Note that these /zq values 
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Figure 2.18: Contributions of individual channels to the total O(af) 6-quark 
production cross section at fixed-target pp experiments versus beam momen- 
tum for rrib = 4.75 GeV/c 2 . Plotted are the contributions of the gg (MS 
scheme, short-dashed line) and qq (DIS scheme, long-dashed line) channels, 
and the absolute value of the contributions of the qg (DIS scheme, upper 
dotted line) and qg (DIS scheme, lower dotted line) channels. 
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Figure 2.19: The /i dependence of the resummed cross section for 6-quark 
production at HERA-B with m& = 4.75 GeV/c 2 for the qq channel in the DIS 
scheme and in the MS scheme, and for the gg channel in the MS scheme. 
Plotted are a T q f (lower solid line DIS scheme, middle solid line MS scheme) 
and a v g e g (upper solid line). Also we plot the sum + crW | app | app 
(lower dotted line for qq in the DIS scheme, middle dotted line for qq in the 
MS scheme, and upper dotted line for gg). 
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are not exactly the same as those used in section 2.2, where Ho/m t = 10% and 
A i o/ m t = 25% for the qq and gg channels respectively, which predicted the mass 
dependence of the top quark cross section. The /xq parameters there were again 
chosen via the criterion that the higher order terms in the perturbation theory 
should not be too large. 

It is illuminating to compare fig. 2.19 with the corresponding plot for the 
top quark case fig. 2.11. There one can infer that if we take the slightly larger 
fio values given above there is very little change in the top quark cross section. 
The reason is that in this case the gg channel makes only a small contribution 
and the /j,q dependence in the qq channel reflects the small variation of the 
running coupling constant at a scale fi = 175 GeV/c 2 . As the running coupling 
constant varies more rapidly at a scale /i = 4.75 GeV/c 2 , the /io parameters 
should be taken from measurements at the lower scale and then used in the 
prediction of the top quark cross section. This emphasizes the importance 
of the proposed measurement at HERA-B. It is clear from fig. 2.19 that we 
cannot choose Ho/mb = 25% for the gg channel for bottom quark production 
but we can choose Ho/m t = 36% for the gg channel for top quark production, 
with very little change in the value of the top quark cross section. Both sets 
of parameters yield cross sections which are within the error bars of the recent 
CDF [U| and DO [|14| experimental results for the top quark cross section. 
Therefore our cut off parameters do have experimental justification. We would 
also like to point out that an application of the principal value resummation 
method has been recently completed by Berger and Contopanagos || leading 
to essentially the same mass dependence of the top cross section as reported in 
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||, which again justifies our choice for /i . Finally note that we could just as 
easily have chosen to work in the MS scheme for both channels by changing /xq 
in the qq channel to /io ~ 1.3 GeV/c 2 . The reason the DIS scheme is preferred 
is simply because it has a larger radius of convergence. 

Using the values of /io that we chose from the previous graphs, we proceed 
to plot the resummed cross section for 6-quark production at fixed-target pp 
experiments versus beam momentum. We present the results in fig. 2.20 for 
the qq and gg channels. For comparison the exact NLO results are shown as 
well. The resummed cross sections were calculated with the cut > Sq while 
no such cut was imposed on the NLO result. Then, in fig. 2.21 we plot the 
resummed cross section and the improved total cross section (2.2.6) (where we 
have also taken into account the small negative contributions of the qg and qg 
channels) versus beam momentum and, for comparison, the total exact NLO 
cross section for the three choices // = mt > /2, m&, and 2mj,. The total NLO 
cross section for 6-quark production at HERA-B (beam energy 820 GeV) is 
28.8 nb for \x = m&/2; 9.6 nb for /j = mj; and 4.2 nb for /j = 2m&. The 
resummed cross section is 18 nb. The improved total cross section for 6-quark 
production at HERA-B is 19.4 nb. 

2.4 Conclusions 

We have presented NLO and resummed results for the cross sections for 
top quark production at the Fermilab Tevatron and for bottom quark produc- 
tion at HERA-B and at fixed-target pp experiments in general. In both cases 
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Figure 2.20: Resummed and NLO cross sections versus beam momentum for 
6-quark production at fixed-target pp experiments for m b = 4.75 GeV/c 2 . 
Plotted are the resummed cross sections for the qq channel in the DIS scheme 
for /j, = 0.6 GeV/c 2 (short-dashed line), and for the gg channel in the MS 
scheme for /x = 1.7 GeV/c 2 (long-dashed line); and the 0(a^) cross sections 
for the gg channel in the MS scheme and the qq channel in the DIS scheme 
(upper and lower solid lines, respectively). 
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Figure 2.21: Resummed, improved, and NLO cross sections versus beam mo- 
mentum for 6-quark production at fixed-target pp experiments for m h = 4.75 
GeV/c 2 . Plotted are the total resummed cross section (long-dashed line), the 
total improved cross section (short-dashed line), and the total 0(a^) cross 
section (// = mj solid line, \i = to&/2 upper dotted line, \i = 2mb lower dotted 
line) . 
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we found that the threshold region gives the main contribution to the NLO 
cross sections. Approximations for the soft gluon contributions in that region 
have been compared with the exact results. The resummation of the leading 
S+V logarithms produces an enhancement of the NLO results. 

For top quark production at the Fermilab Tevatron we saw that the qq 
channel is dominant. We found that the resummation for this channel is 
relatively well behaved and that the scale dependence of the NLO cross section 
is relatively flat. The total NLO cross section for top quark production with 
m t = 175 GeV/c 2 at the Fermilab Tevatron with = 1.8 TeV is 4.8 pb. The 
upper and lower values of the improved cross section for top quark production 
at the Fermilab Tevatron are 5.8 pb and 5.1 pb, respectively. 

For bottom quark production at fixed-target pp experiments it was shown 
that the gg channel is dominant. The leading S+V approximation is not very 
good in the gg channel in the MS scheme in the kinematic region that is 
important for bottom quark production at HERA-B. The addition of sublead- 
ing S+V terms and Coulomb terms improves the approximation considerably. 
The resummation is not as successful as in the qq channel and the scale depen- 
dence of the NLO cross section is much bigger. The total NLO cross section 
for 6-quark production at HERA-B (beam energy 820 GeV) with nib = 4.75 
GeV/c 2 is 9.6 nb (for /j, = m b ). The improved total cross section for 6-quark 
production at HERA-B is 19.4 nb. 



53 



Bibliography 

[1] W. Beenakker, W. L. van Neerven, R. Meng, G. A. Schuler, and J. Smith, 
Nucl. Phys. B351, 507 (1991). 

[2] R. Meng, G. A. Schuler, J. Smith, and W. L. van Neerven, Nucl. Phys. 
B339, 325 (1990). 

[3] W. Beenakker, H. Kuijf, W. L. van Neerven, and J. Smith, Phys. Rev. D 
40, 54 (1989). 

[4] E. Laenen, J. Smith, and W. L. van Neerven, Nucl. Phys. B369, 543 

(1992) . 

[5] A. D. Martin, R. G. Roberts, and W. J. Stirling, Phys. Lett. B 306, 145 

(1993) . 

[6] H. Plothow-Besch, 'PDFLIB: Nucleon, Pion and Photon Parton Density 
Functions and a s Calculations', Users's Manual - Version 4.16, W5051 
PDFLIB, 1994.01.11, CERN-PPE. 

[7] P. Nason, S. Dawson, and R. K. Ellis, Nucl. Phys. B303, 607 (1988). 



54 

[8] E. Laenen, J. Smith, and W. L. van Neerven, Phys. Lett. B 321, 254 
(1994). 

[9] E. Berger and H. Contopanagos, Phys. Lett. B 361, 115 (1995); ANL- 
HEP-PR-95-82, hep-ph/9603326, 1996. 

[10] S. Catani, M. L. Mangano, P. Nason, and L. Trentadue, CERN-TH-96-86, 
hep-ph/9604351, 1996. 

[11] H. Aubrecht et al., An Experiment to Study CP violation in the B system 
Using an Internal Target at the HERA Proton Ring, Letter of Intent, 
DESY-PRC 92/04 (1992). 

[12] H. Aubrecht et al., An Experiment to Study CP violation in the B system 
Using an Internal Target at the HERA Proton Ring, Progress Report, 
DESY-PRC 93/04 (1993). 

[13] (CDF Collaboration) F. Abe et al, Phys. Rev. Lett. 74, 2626 (1995). 

[14] (DO Collaboration) S. Abachi et al, Phys. Rev. Lett. 74, 2632 (1995). 



Chapter 3 

Top and bottom quark inclusive differential 
distributions 

The inclusive transverse momentum and rapidity distributions for top 
quark production at the Fermilab Tevatron and bottom quark production at 
HERA-B are presented both in order in QCD and using the resummation 
of the leading soft gluon corrections in all orders of QCD perturbation theory. 
The resummed results are uniformly larger than the 0(a^) results for both 
distributions. 

3.1 Introduction 

At the Tevatron, the top quark is mainly produced through tt pair pro- 
duction from the light mass quarks and gluons in the colliding proton and 
antiproton. Both the top quark and the top antiquark then decay to (W, b) 
pairs, and each W boson can decay either hadronically or leptonically. The 
6-quark becomes an on-mass-shell 5-hadron which subsequently decays into 
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leptons and (charmed) hadrons. A large effort is being made to reconstruct the 
top quark mass from the measured particles in the decay, which is complicated 
by the fact that the neutrinos are never detected. Also there are additional 
jets so it is not clear which ones to choose to recombine [I], £J . The best chan- 
nel for this mass reconstruction is where both W bosons decay leptonically, 
one to a (e, u e ) pair, the other to a (//, z/ M ) pair (a dilepton event) because the 
backgrounds in this channel are small. When only a single lepton is detected 
then it is necessary to identify the b quark in the decay to remove large back- 
grounds from the production of W+ jets ||. In all cases the reconstruction 
of the particles in the final state involves both the details of the production of 
the top quark-antiquark pair as well as the knowledge of their fragmentation 
and decay products. 

In the analysis of the decay distributions one needs knowledge of the 
inclusive differential distributions of the heavy quarks in transverse momentum 
Pt and rapidity Y. These distributions are known in NLO j|, We present 
an analysis of the resummation effects on the inclusive transverse momentum 
distribution of the top quark assuming that it has a mass of 175 GeV/c 2 . Also 
we discuss here how the resummation effects modify the rapidity distribution 
of the top quark. Since there have been suggestions of using the mass and 
angular distributions in top quark production to look for physics beyond the 
standard model || it is very important to know the normal QCD predictions 
for these quantities. 

We also present a corresponding analysis for the pr and Y distributions 
for bottom quark production at the HERA-B experiment. 
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3.2 Soft gluon approximation to the inclusive 
distributions 

The partonic processes under discussion will be denoted by 

i(h)+j(k 2 ) - Q(pi)+Q(p2) + g(k 3 ), (3.2.1) 

where i,j = g,q,q. The kinematical variables 

s=(k 1 + k 2 ) 2 , t x = (k 2 - p 2 f - m 2 , Ul = (h - p 2 ) 2 - m 2 ,(3.2.2) 

are introduced in the calculation of the corrections to the single particle inclu- 
sive differential distributions of the heavy (anti) quark. We do not distinguish 
in the text between the heavy quark and heavy antiquark since the distribu- 
tions are essentially identical in our calculations. Here s is the square of the 
parton-parton cm. energy and the heavy quark transverse momentum is given 
by Pt = (tiUi/s — m 2 ) 1 / 2 . The rapidity variable is defined by exp(2y) = u\jt\. 
Also as before we define s 4 = s + t± + ui. 

The transverse momentum px of the heavy quark is related to our previous 
variables by 

h = -l{ s -s 4 - [(s - s A ) 2 - Asm 2 ?} 1 ' 2 } , (3.2.3) 
«i = ~{s -s 4 + [(s - s A ) 2 - Asm 2 ,] 1 / 2 } , (3.2.4) 
with my = m 2 + p 2 T . The double differential cross section is therefore 



2 d 2 a ij (s,t 1 ,u 1 ) _ 2 , 2] i/ 2 d 2 (Jij(s, s 4 ,Pt) /„ 9fl 

8 dh d Ul ~ S[{S ~ S4) " 4Smd dp 2 T ds 4 ' {6 - 2 - b) 
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with the boundaries 



<Pt < --m 2 , < s 4 < s - 2m T s 1/2 . 



(3.2.6) 



The 0(a k s ) contribution to the inclusive transverse momentum distribution 
ddij/dpT is given by 

d<rl-\s,p 2 T ) 2 



2fc-l 



rs-2m T s 1 / 2 



dp\ 



s ;=0 JO 



x 



X 



1 



B { 2 \ 



(3.2.7) 



where we have inserted an extra factor of 2 so that / dp\ da / ' dp\ = a tot . After 
some algebra we can rewrite this result as 



da j (S d20 2fe_1 r /■s-2m T s 1 /2 



ds4 — In' — ^ 

s 4 m 2 



x 



dffi? (s, s 4 , Pt) dalj' (s, 0, 



(0), 



dp\ 



dp\ 



} 



/ + 1 V m 2 > dpi J ' V ; 



with the definition 

dafMs, s 4 ,p T ) 



dp 2 - 



ro-J(s,s 4 ,Pr) , 



(3.2.9) 



s[(s - s 4 ) 2 - Asm 2 ,} 1 / 2 lJ 
where daf^ l (s,0,p T ) / 'dp\ = da^\s,p T )/dp T again represents the Born differ- 
ential pt distribution. For the qq and gg subprocesses we have the explicit 
results 

1 



da qq (s,s 4 ,p 2 T ) = 2v ^) Kq9 NC F - 



dp\ 



s [(s - s 4 ) 2 - Asm 2 ,] 1 / 2 



x 



(s — s 4 ) 2 — 2sp: 



T 



(3.2.10) 
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and 



daf g \s,s Al p 2 T ) 
dp? 



A7rai(n 2 )K gg NC F 



s [(s - s 4 ) 2 - 4sm|]V 2 



x 



x 



C F -C. 



Tfln 



(s ~ 54)2 ; 2sm ^ + (3.2.11) 



Since the above formulas are symmetric under the interchange ii <-> Mi the 
heavy quark and heavy antiquark inclusive pr distributions are identical. Note 
that (3.2.8) is basically the integral of a plus distribution together with a 
surface term. 

The corresponding formula to (3.2.8) for the rapidity y of the heavy quark 
is obtained by using 

*i = - ~ S4) (l-tanhy), (3.2.12) 



Ul = _if_£l2(i + tanhy) 



(3.2.13) 



The double differential cross section is therefore 



2 d 2 a ij (s,t 1 ,u 1 ) 



dt\ du\ 



2s 



2 cosh 2 y d 2 a ij (s, s 4 ,y) 
s — S4 dy ds4 



(3.2.14) 



with the boundaries 



"^(l^f) <y< ^ ln (l^|) ' < s 4 < s -2ms 1 / 2 cosh (3.2.15) 

where f3 2 = 1 — 4m 2 /s. The 0(a k s ) contribution to the inclusive rapidity 
distribution daij/dy is given by 



dy 



2fc-l 



-■- M-2ms'/ 2 cosh j/ 
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x{-W (4)6>(* - A) + -1- (^)S(s 4 )} 

L s 4 ^m^ 7 / + 1 ^m 1 ' > 

After some algebra we can rewrite this result as 

d(T^(s V) 2fe_1 r /•s-2ms 1 /2 coshy i 



dy 



x 



1=0 



with the definition 



1 

dy dy 

1 ln m f s - 2mgl/2 cosh V \ d °j? ( s > 1 (3217) 
^ m 2 ' dy ^ 



l + l 



dalf(s,s 4 ,y) _ 



s — S4 
2s 2 cosh 2 ?/ 



(3.2.18) 



dy 

where dalj\s,0,y)/dy = da^\s,y)/dy again represents the Born differential 
1/ distribution. For the qq and yy subprocesses we have the explicit formulas 

da$(s,s A ,y) 2 2 s - s 4 
= ™> )K q ,NC F2s2cQsh2y 



x - (S 2 ^ (cosh 2 y + sinh 2 y) + ^] , (3.2.19) 
L 2s 2 cosh y v 7 s J 



and 



dy '99 * 2s 2 cosh 2 ?/ 

x ^-^77^2- x c 



— ^A—r. 

4s z cosh y 
8m 2 s cosh 2 y / Am 2 s cosh 2 y \ 



cosh 2 y + sinh 2 y 



- s 4 ) 2 



-s 4 ) 2 



(3.2.20) 
the 



Since the above formulas are symmetric under the interchange t 1 <-> 
heavy quark and heavy antiquark inclusive y distributions are identical. Also 
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(3.2.17) is again of the form of a plus distribution together with a surface term. 
Finally, we note that the terms in (3.2.8) and (3.2.17) are all finite. 

3.3 Resummation procedure in parton-parton 
collisions 



We now consider the resummation of the order a s contributions to the pr 
distribution. We have 

daij(s,p\) y, dalj\s,p^) 
dp 2 T ^ dp 2 T 

s -2m T s 1 / 2 df(s 4 /m 2 ,m 2 /fj J 2 ) 

uS 4 — 

so US4 



I ■ dermis, s 4 ,pl) _ dermis, 0,pp 
>- dpi dpi 

<s-2m T s 1 / 2 m 2 x da§\s, p 2 T ) 



+i \-^ 2 — v) dpi • (3 - 3 - 1} 

Note that, as in chapter 2, we now have cut off the lower limit of the S4 
integration at S4 = sq because a s diverges as S4 — > 0. This parameter so must 
satisfy the conditions < s < s — 2m T s 1 ^ 2 and s /m 2 « 1. The derivative 
of f(s4/m 2 ,m 2 /fi 2 ) is obtained from (2.1.24). It is equal to 

df(si/m 2 ,m 2 / ji 2 ) 



l{2^a s (^,m 2 )ln^ + r ? } 



d$4 S4 ^ 7r m 2 m 2 

x exp {A^a s 4, m 2 ) In 2 -L} ^ 4/ J 
L 7T m z m 2 > r(l + 77) 

x exp(-r ?7£ ) , (3.3.2) 



where we have neglected terms which are higher order in a s . 
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The analogous formula for the rapidity distribution is 

d<Tij(s,y) _ ™da§\s,y) 



dy fo d v 

-L 



s^ms 1 / 2 coshy df (s 4 / 171 2 , 171 2 / jJ, 2 ) 



sq ds4 
f dalf{s,S4,y) _ dalfjs, 0, y) ^ 
dy dy ' 

.(S- 2ms 1 / 2 cosh y m 2 n da$ (s , y) 

+/( 2 '—J — a ' ( 3 - 3 - 3 ) 

v m z fi z ' dy 

with the conditions < So < s — 2ms 1 / 2 coshy and So/m 2 « 1. 

3.4 Top quark differential distributions 

Since the px distribution in hadron-hadron collisions is not altered by 
the Lorentz transformation along the collision axis from the parton-parton 
cm. frame, we can write an analogous formula to (2.2.1) for the heavy-quark 
inclusive differential distribution in p\ 

d^(S, m 2 ,p 2 T ) ,i da§\TS,m 2 ,p 2 T ^ 2 ) 

dpi '~ 4r-^4>- ''in '~ ' 

In the case of the all-order resummed expression the lower boundary in (3.4.1) 
has to be modified according to the condition s < s — 2m T s 1 / 2 (see section 
3.3). The all-order resummed differential distribution in p\ is 

daW(S,m 2 ,p 2 T ) x . f 1 ^ da tJ {rS,m 2 ,p 2 T ,^ 2 ) 

— dpj — = dT ^ {r ^ } — m — ' ( } 

with daij/dpl given in (3.3.1) and 

(m T + (m 2 T + SQ) 1 / 2 ) 2 
t = . (3.4.3) 
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The corresponding formula to (3.4.1) for the heavy quark inclusive differential 
distribution in Y is 

<>(W.y) H (Tj(lI) ^«rf. (344) 

dY *-r? Jim 2 cosh 2 y/S dy 

Order by order in perturbation theory the heavy quark rapidity plots in the 
parton-parton cm. frame show peaks away from y = (see fig. 7 in 0). 
However, upon folding with the partonic densities the heavy quark rapidities 
in the hadron-hadron cm. frame peak near Y = 0. Therefore we will assume 
that the plots for the resummed rapidity distribution show a similar feature. 
The all-order resummed differential distribution in Y is therefore taken to be 

doW(S,m* t Y) = ji dr ^^ d^rS^yjQ ^ 

with doijjdy given in (3.3.3) and 

(m cosh y+ (m 2 cosh 2 y + sq) l / 2 ) 2 
t = . (3.4.6) 

The hadronic heavy quark rapidity Y is related to the partonic heavy quark 
rapidity y by 

Y = y + \\n^. (3.4.7) 

We now specialize to top quark production at the Fermilab Tevatron where 
\/~S = 1.8 TeV and choose the top quark mass to be m t = 175 GeV/c 2 . 
In the presentation of our results for the exact, approximate and resummed 
hadronic cross sections we use the same (MRSD_' ) parametrization for the 
parton distributions as we did in the previous chapter. 
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Since we know the exact 0(a^) result, we can make an even better es- 
timate of the differential distributions by calculating the perturbation theory 
improved pt and Y distributions. We define the improved pr distribution by 

drr imp da TCS drr W drr W 

aa H _ aa H _,_ aa H I aa H I fo a o\ 

i t ' i exact , app j yo.Q.O) 

dpx dpr dpT dpx 
and the improved Y distribution by 

d^_d_aW ,do$, dag , 

dY ~ dY + dY |cxact dY |app ' 1 ' 

to exploit the fact that da^/dpr | exact and da^ /dY | cxact are known and 
d<jjp/dpT |ap P and daff/dY | app are included in da r ff/dpT and da T ff/dY re- 
spectively. We note that here da^ denotes the 0(a™ +2 ) contribution to the 
differential cross section. Moreover , da^> | exact denotes the exact calculated 
differential cross section, and da^ | app the approximate one where only the 
leading soft gluon corrections are taken into account. 

First we present the differential pr distributions at =1.8 TeV for 
a top quark mass m t = 175 GeV/c 2 . For these plots the mass factorization 
scale is not everywhere equal to m t . We chose p, — m t in so, fk(s^/m 2 , m 2 //x 2 ) 
and a s , but p = m T in the MRSD_' parton distribution functions and the 
running coupling constant a s (p 2 ). It should be noted, however, that it makes 
little difference if we choose /i — m t everywhere. The difference in the cross 
section is only a few percent so that the changes due to scale dependence are 
insignificant compared with the changes due to higher order resummation. We 
begin with the results for the qq channel in the DIS scheme. In fig. 3.1 we 
show the Born term daff /dpr, the first order exact result da$ / dpx (exact, 
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Figure 3.1: The top quark p T distributions daff jdpT for the qq channel in the 
DIS scheme for a top quark mass m t = 175 GeV/c 2 . Plotted are daf jdp T 
(upper solid line), da^/dpr | exact (lower solid line), da^/dpr | apP (upper 
dotted line), daff/dpT | app (lower dotted line), and da^/dpx (a*o — 0.05 m t 
upper dashed line and /xq = 0.1 m t lower dashed line). 
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the first order approximation da H /dpr | apP , the second order approximation 
da H /dpr | app , and the resummed result da^/dpr for /i = 0.05 m t and for 
fj, — 0.1 m t . These are the same values for /i that we used in chapter 2. As 
we decrease /x the differential cross sections increase. 

We continue with the results for the gg channel in the MS scheme. The 
corresponding plot is given in figure 3.2. In this case the values of /xo have been 
chosen to be /iq = 0.2m t and /io = 0.25m t , again as in chapter 2. The first and 
second order corrections in the gg channel in the MS scheme are larger than 
the respective ones in the qq channel in the DIS scheme. In fact, for the range 
of pt values shown the second-order approximate correction is larger than the 
first-order approximation. Hence, the relative difference in magnitude between 
the improved da^/dpr and the exact O(a^) results is significantly larger than 
that for the qq channel in the DIS scheme. 

We finish our discussion of the differential pt distributions with the results 
of adding the qq and gg channels. The plot appears in figure 3.3. We also 
show the total improved and 0(a^) distributions in fig. 3.4. It is evident that 
resummation produces an enhancement of the exact O(a^) result, with very 
little change in shape. 

Now we turn to a discussion of the differential Y distributions at = 
1.8 TeV for a top quark mass m t = 175 GeV/c 2 . In this case we set the 
factorization mass scale equal to m t everywhere. We begin with the re- 
sults for the qq channel in the DIS scheme. In fig. 3.5 we show the Born 
term daff/dY, the first order exact result da^/dY | CX act, the first order 
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Figure 3.2: The top quark p T distributions da H /dp T for the gg channel in 
the MS scheme for a top quark mass m t = 175 GeV/c 2 . Plotted are daP/dp T 
(upper solid line at large pt), da^ /dpr | exact (lower solid line at large pt), 
da^/dpT | app (lower dotted line), daff/dpT | apP (upper dotted line), and 
da r ff jdpT (/io = 0.2 mt upper dashed line and /xq = 0.25 m t lower dashed line). 
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Figure 3.3: The top quark p T distributions da^ /dp T for the sum of the qq and 
gg channels for a top quark mass m t = 175 GeV/c 2 . Plotted are /dpr 
(upper solid line), daff/dpr | exact (lower solid line), da^/dpr | apP (upper 
dotted line), dajp /dpx | app (lower dotted line), and do^jdpT (upper and 
lower dashed lines). 
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Figure 3.4: The top quark p T distributions d<r H /dpT for the sum of the 
qq and gg channels for a top quark mass m t = 175 GeV/c 2 . Plotted are 
d(jjp /dpT + da^ /dpT | exact (solid line) and da^ v /dpx (upper and lower dashed 
lines) . 
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Figure 3.5: The top quark Y distributions da H jdY for the qq channel in the 
DIS scheme for a top quark mass m t = 175 GeV/c 2 . Plotted are daf jdY 
(upper solid line), da^ jdY | cxact (lower solid line), da^ jdY | app (upper dotted 

(2) 

line), da H jdY | app (lower dotted line), and da^/dY (/i = 0.05 m t upper 
dashed line and /io = 0.1 m t lower dashed line). 
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approximation da H jdY | app , the second order approximation da H jdY | app , 
and the resummed result da^/dY for fio = 0.05 m t and fio = 0.1 m t . 

We continue with the results for the gg channel in the MS scheme. The 
corresponding plot is given in figure 3.6. Here, the values of /z are fj, = 0.2 m t 
and fj,o = 0.25 m t . As in the case of the px distributions, the first and second 
order corrections in this channel are larger than the respective ones in the qq 
channel in the DIS scheme. For the range of Y values shown the second-order 
approximate correction is larger than the first-order approximation. Again, 
as in the pt distributions, the relative difference in magnitude between the 
improved da^/dY and the exact 0(a 3 ) results is significantly larger than 
that in the qq channel in the DIS scheme. 

Finally, we conclude our discussion of the differential Y distributions by 
showing the results of adding the qq and gg channels. The plots appear in 
figures 3.7 and 3.8. Again, it is evident that resummation produces a non- 
negligible modification of the exact 0(a 3 s ) result. However, the shape of the 
distribution is unchanged. 



3.5 Bottom quark differential distributions 

In this section we present some results on the inclusive transverse mo- 
mentum and rapidity distributions of the bottom quark at HERA-B. 

We begin with the px distributions. For these plots the mass factorization 
scale is not everywhere equal to m b . We chose \i = m b in s , /fc(s 4 /m 2 , m 2 //x 2 ) 
and a s , but jx = mr in the MRSD_' parton distribution functions and the 
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Figure 3.6: The top quark Y distributions da$/dY for the gg channel in the 
MS scheme for a top quark mass m t = 175 GeV/c 2 . Plotted are da^/dY 
(upper solid line at Y = 0), da^/dY | cxact (lower solid line at Y = 0), 
da^/dY | app (lower dotted line), daff/dY | app (upper dotted line), and 
d(jff S /dY (/x = 0.2 m t upper dashed line and /iq = 0.25 m t lower dashed 
line) . 
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Figure 3.7: The top quark Y distributions da H jdY for the sum of the qq 
and gg channels for a top quark mass m t = 175 GeV/c 2 . Plotted are daf jdY 
(upper solid line), da^ jdY | cxact (lower solid line), da^ jdY | app (upper dotted 
line), duff jdY | app (lower dotted line), and da^/dY (upper and lower dashed 
lines) . 
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Figure 3.8: The top quark Y distributions dan/dY for the sum of the qq 
and gg channels for a top quark mass m t = 175 GeV/c 2 . Plotted are 
dajp I 'dY + da^ / 'dY | cxact (solid line) and da^ v /dY (upper and lower dashed 
lines). 
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running coupling constant a s (/i 2 ). In fig. 3.9, we give the results for the qq 
channel in the DIS scheme. We plot the Born term daff/dpr, the first order 
exact result da$ /dp T \ exact , the first order approximation da^/dpr \ ap p, the 
second order approximation daff /dpx | apP , and the resummed result da^ s jdpT 
for [i = 0.6 GeV/c 2 . This is the same value for /x that we used in chapter 
2 for the total cross section. If we decrease the differential cross sections 
will increase. The resummed distribution was calculated with the cut s 4 > s 
while no such cut was imposed on the phase space for the individual terms in 
the perturbation series. We continue with the results for the gg channel in the 
MS scheme. The corresponding plot is given in fig. 3.10. We note that the 
corrections in this channel are large. In fact the exact first-order correction 
is larger than the Born term and the approximate second-order correction is 
larger than the approximate first-order correction. In this case the value of 
Hq has been chosen to be fj, = 1.7 GeV/c 2 as in chapter 2. In fig. 3.11 
we plot the improved px distribution for the sum of all channels, where we 
have included the small negative contributions of the qg and qg channels. For 
comparison we also show the total exact NLO results for \i = rrib/2, nib, and 
2m b . The improved px distribution is uniformly above the exact 0(a^) results. 
We see that the effect of the resummation exceeds the uncertainty due to scale 
dependence. 

We finish with a discussion of the Y distributions. In this case we set 
the factorization mass scale equal to mj everywhere. We begin with the qq 
channel. In fig. 3.12 we show the Born term da^/dY, the first order exact 
result d(jft jdY | cxac t, the first order approximation da^ jdY | app , the second 
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Figure 3.9: The bottom quark p T distributions da H /dpx at HERA-B for the 
qq channel in the DIS scheme for m b = 4.75 GeV/c 2 . Plotted are daff /dpr 
(upper solid line), da^/dpr | CX act (lower solid line), da^/dpr | apP (upper 
dotted line), da^/dpT | app (lower dotted line), and da^/dpT for /i = 0.6 
GeV/c 2 (dashed line). 
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Figure 3.10: The bottom quark p T distributions da$ jdpT at HERA-B for the 
gg channel in the MS scheme for m b = 4.75 GeV/c 2 . Plotted are /dpr 
(lower solid line), da$ /dpr | exact (upper solid line), da$ /dpr | apP (lower dot- 
ted line), daff /dpT | apP (upper dotted line), and da r ^ s /dpx for p, = 1.7 GeV/c 2 
(dashed line). 
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Figure 3.11: The bottom quark p T distributions da H /dp T at HERA-B 
for the sum of all channels for mt = 4.75 GeV/c 2 . Plotted are 
dajp /dpT + dcrff* jdpT | C xact (/U = m b solid line, n = m b /2 upper dotted line, 
ji = 2rrib lower dotted line) and da^ v /dpr (dashed line). 
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Figure 3.12: The bottom quark Y distributions da H jdY at HERA-B for the 
qq channel in the DIS scheme for m h = 4.75 GeV/c 2 . Plotted are do^ jdY 
(upper solid line), da^ jdY | exact (lower solid line), da^ jdY | app (upper dot- 
ted line at positive Y), do^ jdY | app (lower dotted line at positive Y), and 
da^/dY for ^ = 0.6 GeV/c 2 (dashed line). 
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order approximation da H jdY | app , and the resummed result da^/dY for 
Ho = 0.6 GeV/c 2 . Again, the resummed distribution was calculated with 
the cut S4 > so while no such cut was imposed on the phase space for the 
individual terms in the perturbation series. We continue with the results for 
the gg channel in the MS scheme. The corresponding plot is given in fig. 3.13. 
Here, the value of /i is /io = 1-7 GeV/c 2 . The corrections in this channel 
are large as was the case for the px distributions. In fig. 3.14 we plot the 
improved Y distribution for the sum of all channels, where we have included 
the small negative contributions of the qg and qg channels. For comparison 
we also show the total exact NLO results for /j, = m b /2, m b , and 2m b . The 
improved Y distribution is uniformly above the exact 0(a^) results. Again, 
we see that the effect of the resummation exceeds the uncertainty due to scale 
dependence. 

3.6 Conclusions 

We have shown that the resummation of soft gluon radiation produces a 
small difference between the perturbation theory improved distributions in p T 
and Y and the exact 0(a%) distributions in px and Y for the qq reaction in 
the DIS scheme for the values of fi chosen. However, for the gg channel in 
the MS scheme the resummation produces a large difference. The difference 
between the resummed and the exact 0(a^) distributions depends on the mass 
factorization scheme (DIS or MS), the factorization scale //, as well as the 
specific reaction under consideration (qq or gg). For top quark production 
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Figure 3.13: The bottom quark Y distributions da H jdY at HERA-B for the 
gg channel in the MS scheme for m b = 4.75 GeV/c 2 . Plotted are do^ jdY 
(lower solid line), da^ jdY | cxact (upper solid line), da^ jdY | app (lower dotted 
line), duff /dY | app (upper dotted line), and da^/dY for /i = 1.7 GeV/c 2 
(dashed line). 
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Figure 3.14: The bottom quark Y distributions da H /dY at HERA-B 
for the sum of all channels for mt = 4.75 GeV/c 2 . Plotted are 
da^x ' /dY + da^H jdY | cxact (/i = m b solid line, \x = m b /2 upper dotted line, 
/i = 2m& lower dotted line) and da^ p /dY( dashed line). 
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at the Fermilab Tevatron with m t = 175 GeV/c 2 the gg channel is not as 
important numerically as the qq channel. However, since the corrections for the 
gg channel are quite large, resummation produces a non-negligible difference 
between the perturbation theory improved and the exact 0(a^) distributions 
when adding the two channels. However, the shapes of the distributions are 
essentially unchanged. For bottom quark production at HERA-B with m& = 
4.75 GeV/c 2 the gg channel is dominant so the enhancement of the NLO 
distributions is much bigger. 
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Chapter 4 

Resummation of singular distributions in 
QCD hard scattering 



We discuss the resummation of distributions that are singular at the elas- 
tic limit of partonic phase space (partonic threshold) in QCD hard-scattering 
cross sections, such as heavy quark production. We show how nonleading 
soft logarithms exponentiate in a manner that depends on the color structure 
within the underlying hard scattering. This result generalizes the resumma- 
tion of threshold singularities for the Drell-Yan process, in which the hard 
scattering proceeds through color-singlet annihilation. We illustrate our re- 
sults for the case of heavy quark production by light quark annihilation and 
gluon fusion, and also for light quark production through gluon fusion. 



4.1 General formalism 



In hard scattering cross sections factorized according to perturbative QCD 
the calculable short- distance function includes distributions that are singular 
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when the total invariant mass of the partons reaches the minimal value neces- 
sary to produce the observed final state. Such singular distributions can give 
substantial QCD corrections to any order in a s . 

Expressions that resum these distributions to the short-distance func- 
tion of Drell-Yan cross sections to arbitrary logarithmic accuracy have been 
known for some time [I]. It has also been observed that leading distributions, 
and hence leading logarithms in moment space, are the same for many hard 
QCD cross sections. We used this fact as the basis for our estimates of heavy 
quark production cross sections and differential distributions in the previous 
chapters. In chapter 2 we pointed out the inadequacy of the leading log ap- 
proximation, particularly for the gg channel. In this chapter, we shall extend 
our analysis to the level of next-to-leading logarithms. We shall exhibit a 
method by which nonleading distributions may be treated, and will illustrate 
this method in the cases of heavy quark production through light quark anni- 
hilation and gluon fusion, and light quark production through gluon fusion. 

We consider the inclusive cross section for the production of one or more 
particles, with total invariant mass Q. Examples include states produced 
by QCD, such as heavy quark pairs or high-py jets, in addition to massive 
electroweak vector bosons, virtual or real, as in the Drell-Yan process. 

To be specific, we shall discuss the summation of ("plus") distributions, 
which are singular for z = 1, where 




for the production of a heavy quark pair of total invariant mass Q, with s 
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the invariant mass squared of the incoming partons that initiate the hard 
scattering. We shall refer to z = 1 as "partonic threshold" I, or more accurately 
the "elastic limit." We assume that the cross section is defined so that there 
are no uncancelled collinear divergences in the final state. 

The main complications relative to Drell-Yan []J involve the exchange 
of color in the hard scattering, and the presence of final-state interactions. 
In fact, these effects only modify partonic threshold singularities at next-to- 
leading logarithm, and we give below explicit exponentiated moment-space 
expressions which take them into account at this level. At the next level 
of accuracy, we shall see that resummation requires ordered exponentials, in 
terms of calculable anomalous dimensions. 

The properties of QCD that make this organization possible are the fac- 
torization of soft gluons from high-energy partons in perturbation theory ||, 
and the exponentiation of soft gluon effects 0. Factorization is represented 
by fig. 4.1 for the annihilation of a light quark pair to form a pair of heavy 
quarks. In this figure, momentum configurations that contribute singular be- 
havior near partonic threshold are shown in a cut diagram notation [0. As 
shown, it is possible to factorize soft gluons from the "jets" of virtual and 
real particles that are on-shell and parallel to the incoming, energetic light 
quarks, as well as from the outgoing heavy quarks. Soft-gluon factorization 

1 We emphasize that by partonic threshold, we refer to cm. total energy of the in- 
coming partons for a fixed final state; heavy quarks, for example, are not necessarily 
produced at rest. 
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from incoming light-like partons is a result of relativistic limit [[|, while fac- 
torization from heavy quarks, even when they are nonrelativistic is familiar 
from heavy-quark effective theory. Once soft gluons are factored from them, 
the jets may be identified with the parton distributions of the initial state 
hadrons. The hard interactions, labelled Hj and Hj in fig 4.1, corresponding 
to contributions from the amplitude and its complex conjugate, respectively, 
are labelled by the overall color exchange in each. A general argument of how 
the exponentiation of Sudakov logarithms follows from the factorization of soft 
and hard parts and jets is given in ||. 

For example, with the quark-antiquark process shown, the choice of color 
structure is simple, and may, for instance, be chosen as singlet or octet. To 
make these choices explicit, we label the colors of the incoming pair % and j 
for the quark and antiquark respectively, and of the outgoing (massive) pair A; 
and I for the quark and antiquark. The hard scattering is then of the generic 
form 

Hi = />i(Q 2 //A a.G" 2 )) SjiSik , (4.1.2) 

for singlet structure (annihilation of color) in the s-channel. For the s-channel 
octet, or more generally adjoint in color SU(N), we have, analogously 

N 2 -l 

H A = h A (Q 2 /^ 2 ,a s ^ 2 )) £ [Tf) 

c=l 

with TP the generators in the fundamental representation. The functions 
hj are, as indicated, infrared safe, that is, free of both collinear and infrared 
divergences, even at partonic threshold. 

Taking into account possible choices of Hi and H}, an expression that 



T (F) 



Ik 



(4.1.3) 
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Figure 4.1: Cut diagram illustrating momentum configurations that give rise to 
threshold enhancements in heavy quark production: (a) General factorization 
theorem. Away from partonic threshold all singularities in the "short-distance" 
subdiagram H/S cancel; (b) Expanded view of H/S near threshold, showing 
the factorization of soft gluons onto eikonal (Wilson) lines from incoming and 
outgoing partons in the hard subprocess. Hj and Hj represent the remaining, 
truly short-distance, hard scattering. 
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organizes all singular distributions for heavy quark production is 

da hlh2 r dx a dx h 



12 12 [ Xa — <l>a/hA X a,Q 2 )<l>b/h 2 {Xb,Q 2 ) 



dQ 2 d cos 9 dy ~£ Yj J x a Xb 

(4.1.4) 

where y is the pair rapidity and 9 is the scattering angle in the pair center of 



mass frame. The indices / and J label color tensors, such as the singlet ( |4.1.2| ) 
and octet (|4.1.3|) , with which we contract the color indices of the incoming 
and outgoing partons that participate in the hard scattering. The variable S 
is the invariant mass squared of the incoming hadrons. The functions <p a /h are 
parton densities, evaluated at scale Q 2 . The function Q contains all singular 
behavior in the threshold limit, z — *■ 1. Q depends on the scheme in which 
we perform factorization, the usual choices being MS and DIS. Note that the 
resummation may be carried out at fixed y, so long as y is not close to the 
edge of phase space p. 

The color structure of the hard scattering influences contributions to non- 
leading infrared behavior. Not all soft gluons, however, are sensitive to the 
color structure of the hard scattering. Gluons that are both soft and collinear 
to the incoming partons, factorize into the parton distributions of the incom- 
ing hadrons. It is at the level of nonleading purely soft gluons with central 
rapidities that color dependence appears, in the resummation of soft gluon 
effects. Each choice of color structure has, as a result, its own exponentiation 
for soft gluons II. Then, to next-to-leading-logarithm (NLL) it is possible to 
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pick a color basis in which moments with respect to z exponentiate, 

n^(n,y,9,Q 2 ) = f 'dz^nSP&yAa.iQ?)) 

Jo 

= HV J \y,e,Q*)e E "^ , (4.1.5) 
where the color-dependent exponents are given by 



E u (n,e,Q 2 ) = - - {z n - l -l) 

Jo 1 — z 



dlJ g[ ab) [a s ((l-y)(l-z)Q 2 )] 



o 1 -y 



+ g { 2 ab) [a s ((l - z)Q 2 )} +gi I) [a s ((l - z) 2 Q 2 ),6] 
+ g { 3 J >[a s ((l-z) 2 Q 2 ),9] 



(4.1.6) 



The Qi are finite functions of their arguments. The 11^^ are infrared safe 
expansions in a s (Q 2 ). g[ ab ^ and are universal among hard cross sections 
and color structures for given incoming partons a and b, but depend on whether 
these partons are quarks or gluons. On the other hand, #3 summarizes soft 
logarithms that depend directly on color exchange in the hard scattering, and 
hence also on the identities and relative directions of the colliding partons 
(through 8), both incoming and outgoing. 

Just as in the case of Drell-Yan, to reach the accuracy of NLL in the 
exponents, we need g\ only to two loops, with leading logarithms coming 
entirely from its one-loop approximation, and g% and g^p only to a single loop. 
More explicitly, we take 

9 r> = (c„ + a)^ + M^1, (4.1.7) 
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with Ci = Cp (C a) for an incoming quark (gluon), and with K given by 



where rif is the number of quark flavors, g?, is given for quarks in the DIS 
scheme by 



and it vanishes in the MS scheme. As pointed out in |§, one- loop contribu- 
tions to #3 may always be absorbed into the one-loop contribution to #2 and 
the two- loop contribution to g\. Because g% depends upon J, however, it is 
advantageous to keep this nonfactoring process-dependence separate. We shall 
describe how it is determined below. 

First, let us sketch how these results come about J3J. After the normal 
factorization of parton distributions, soft gluons cancel in inclusive hard scat- 
tering cross sections. When restrictions are placed on soft gluon emission, 
however, finite logarithmic enhancements remain, and it is useful to separate 
soft partons from the hard scattering (which is then constrained to be fully vir- 
tual). Soft gluons may be factored from the hard scattering into a set of Wilson 
lines, or ordered exponentials, from which collinear singularities in the initial 
state are eliminated, either by explicit subtractions or by a suitable choice of 
gauge ||. Assuming that the lowest-order process is two-to-two, there will 
be two incoming and two outgoing Wilson lines.i The result, illustrated in 

2 In Drell-Yan and other electroweak annihilation processes, there is a pair of 
incoming lines only. 




(4.1.8) 




(4.1.9) 
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fig. 4.1b, is of the form, H^Sjj, summed over the same color basis as in 
(4.1.4) above. 

The resulting hard scattering and soft-gluon functions both require renor- 
malization, which is organized by going to a basis in the space of color ex- 
changes between the Wilson lines. The renormalization is carried out by a 
counterterm matrix in this space of color structures. For incoming and outgo- 
ing lines of equal masses, such analyses have been carried out to one loop in ||, 
||, and fllPl , and to two loops in a related process in [jllj. For an underlying 



partonic process a + b — > c+d, we then construct an anomalous dimension ma- 
trix T^ b {] cd \ where the indices I and J vary over the various color exchanges 
possible in the partonic process. We write for the renormalization of S 

Sf) = l -[Z s>ir 5jj, + Z s ,jj,6 n ,]S rj ,, (4.1.10) 

where denotes the unrenormalized quantity. The soft function SVj then 
satisfies the renormalization group equation || 

In a minimal subtraction scheme with e = n — 4 

r s (g) = -f ^-Res e _^s(<7,e). (4.1.12) 

This resummation of soft logarithms is analogous to singlet evolution in deeply 
inelastic scattering, which involves the mixing of operators, and hence of par- 
ton distributions. The general solution, even in moment space, is given in 
terms of ordered exponentials which, however, may be diagonalized at leading 
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logarithm. For the resummation of soft logarithms in QCD cross sections, the 
same general pattern holds, with mixing between hard color tensors. Leading 
soft logarithms, however, are next-to-leading overall in moment space, which 
allows the exponentiation (4.1.5) at this level. 



Of course, the analysis is simplest for external quarks, and most compli- 
cated for external gluons. It is also possible to imagine a similar analysis when 
there are more than two partons in the final state. This would be necesary if 
we were to treat threshold corrections to pp — > QQ + jet, for instance, but we 
have not attempted to explore such processes in detail. 

Given a choice of incoming and outgoing partons, next-to-leading loga- 



rithms in the moment variable n exponentiate as in ( [4.1. 5|) in the color tensor 



basis that diagonalizes I^j^ , with eigenvalues A/. The resulting soft func- 
tion g^p is then simply 

#[a s ((l - z) 2 Q 2 ), 9] = -A z [a.((l - z) 2 Q 2 ), 9] , (4.1.13) 

where the eigenvalues are complex in general, and depend on the directions of 
the incoming and outgoing partons, as shown. 

4.2 Applications to qq — » QQ 

These considerations may be illustrated by heavy quark production through 
light quark annihilation, 

q(Pa) + q(Pb) - Q{Pi) + Q{Vi) ■ (4-2.1) 
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In this case, as in elastic scattering g, |TT|, the anomalous dimension matrix 
is only two-dimensional. 

As before we define the invariants 

s = (p a + Pb) 2 , h = {p a ~ Pi) 2 - m 2 , ui = (Pb - Pi) 2 ~ m 2 , (4.2.2) 

with m the heavy quark mass, which satisfy 

s + *! + «! = () (4.2.3) 

at partonic threshold. We also define dimensionless vectors by 

p? = Qv? (4.2.4) 

which obey v 2 = for the light incoming quarks and v 2 = m 2 / Q 2 for the 
outgoing heavy quarks. Note that Q satisfies the kinematic relation 

s = 2Q 2 . (4.2.5) 

We now calculate Fs{g)- The UV divergent 0(a s ) contribution to S is the sum 
of the graphs in fig. 4.2. The counterterms for S are the ultraviolet divergent 
coefficients times our basis color tensors: 

St = CxZs^ii + c 2 Z s ,2i, (4.2.6) 

S 2 = ClZ S ,12 + C 2 Z S ,22- (4.2.7) 

In our calculations we use the axial gauge gluon propagator 

Wlk) = —^N^k), N^(k) = g^- n + , fc +n 2 -^4w, (4.2.8) 
k 2 + le n ■ k (n ■ k) 2 



Figure 4.2: UV divergent one-loop contributions to S for qq 



97 



Figure 4.3: Eikonal rules for qq — > QQ. The gluon momentum flows out of 
the eikonal lines. 

with n M the gauge vector, and eikonal rules for all external lines (fig. 4.3). If 
we denote a typical one-loop correction to cj as u^(SiVi, SjVj, n, Aj, Aj), where 
5 = +1 (—1) for the gluon momenta flowing in the same (opposite) direction 
to the momentum of Vi and A = +1 (—1) for Vi corresponding to a quark 
(antiquark), then we have: 



L 9 J (27T)" 



(27r)™ q 2 + ie { (5iVi ■ q + it) (5jVj ■ q + it) 

AiVi • n P AjVj ■ n P 9 P 

" ■' 1 n 



(SiVi ■ q + ie) (n- q) (SjVj ■ q + ie) (n ■ q) (n ■ q) 2 j ' 

(4.2.9) 

where is a color tensor. P stands for principal value, 



P 11 ! + (-1)' , : * J • (4-2.10) 



(q ■ n) 13 2 \(q ■ n + ie) 13 (—q-n + ie) 13 / 
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We rewrite ( |4.2.9| ) as 

,(-0 



to 



1 l 

h(SiVi, SjVj) - -I 2 (SiVi, n) - -I 2 (5iVi, -n) 



n) - -hiSjVj, -n) + I A (n 2 



(4.2.11) 



where the overall sign is given by 



Sij = Aj Aj 5i Sj. 



(4.2.12) 



We now evaluate the ultraviolet poles of the integrals. For the integrals 
when both Vi and Vj refer to massive quarks we have 



tUV pole 

J l — 


a s 1 T 

L f) > 

7r e 


(4.2.13) 


tUV pole 
L 2 — 


TT e 


(4.2.14) 


T UV pole 

J 3 — 


tt e 


(4.2.15) 


tUV pole 


a s 1 
7T e 


(4.2.16) 



where the Lp is the familiar velocity-dependent eikonal function 



1 - 2m 2 Is / 1-/3 
1 In ■ 



P 



1 + /3 



+ TTl 



(4.2.17) 



with (3 = yl — Am?/ s. The Lj and Lj are complicated functions of the gauge 
vector n. We will see shortly that their contributions are cancelled by the 
inclusion of self energies. Their explicit expressions are: 

1 



L-, - ^[Li(+n) + Lii-n)} 



(4.2.18) 



where 



H+n) 



Vi ■ n 



\J (vi ■ n) 2 — 2m? n 2 j s 
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In 



/ 5(±n) 2m 2 /s 



n\ — \J {vi ■ n) 2 — 2m 2 n 2 / s 



In 



\5(±n) 2m 2 / s — \vi ■ n\ + ^ (vi ■ n) 2 — 2m 2 n 2 / s 
5(±n) n 2 — \vi ■ n\ — J (vi ■ n) 2 — 2m 2 n 2 / s 



(4.2.19) 



5(±n) n 2 — \v{ ■ n\ + <J {vi ■ n) 2 — 2m 2 n 2 /s 
with 5(n) = \v i ■ n\/(vi ■ n). 

When Vi refers to a massive quark and Vj to a massless quark we have 

'vijQ' 



T UV pole 


a s _ 




1 


J l 


2tt ' 




e 


tUV pole 
L 2 


a s 








7T 


e 




tUV pole 


a s ( 




1 


J 3 


2tt ' 


^~ 


e 


T UV pole 




1 




1 4 




■> 





7 + 21n 



t 



1h(4tt) 



[7 + ln(z/j) - ln(47r)] 



7r e 



where 



X ■ ™) 5 



and Vij = Vi ■ Vj. Note that the double poles cancel. 

Finally when both Vi and Vj refer to massless quarks we have 



UV pole 



T UV pole 
L 1 



T UV pole 
J 3 



T UV pole 
J 4 



7T 
Os_ 

2vr 

a s 

2^ 

at* 1 



7 + In [f) - ln(4vr) 
-I[ T + l n (^)-hi(47r)]| 



71 e 



(4.2.20) 

(4.2.21) 
(4.2.22) 
(4.2.23) 

(4.2.24) 



(4.2.25) 
(4.2.26) 
(4.2.27) 
(4.2.28) 



Again, note that the double poles cancel. 

Our calculations are most easily carried out in a color tensor basis con- 
sisting of singlet exchange in the s and u channels, 



Cl = SabS, 



13 ■ 



c 2 — SajSbi 



(4.2.29) 
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The color indices for the incoming quark and antiquark are a and 6, respec- 
tively, and for the outgoing quark and antiquark j and i, respectively. 
In the basis ( |4.2.29|) we find 



a. 



1 5,11 

Ts^i 
r^.22 



c f 



71 



In ( ^jf) ~ L < 3 ~\ \n(v a v b ) - Li - Lj + 2 - iri 



N 



-r 



5,21; 



2vr D 
2tt 

7T 



\VaiVbjJ 

In f - ln(> ai t; bi ) - + In l^- 



1X1 



ln(v aj v bi ) - In 



2m 



- - hi(v a v b ) - Li- Lj + 2 
s y 2 



-^r 5il2 . 



(4.2.30) 



The matrix depends, as expected, on the directions of the Wilson lines, which 
may be reexpressed in terms of ratios of kinematic invariants for the partonic 
scattering. We eliminate the gauge dependence of the heavy quarks by in- 
cluding the self-energy graphs in fig. 4.4. The contribution of the self-energy 
graphs (in the diagonal elements only) is the following: 



71 



(4.2.31) 



Next we absorb the k^k u contribution to #2, since it also appears in Drell Yan. 
This gives a — (a s /7r)CV in the diagonal elements. We also have in the axial 
gauge v a = 1/2. Then in terms of our invariants s, ti, and Ui, the anomalous 
dimension matrix becomes 

C A 



5,11 



0:. s 



7T 



c f 



In 



1ii 



1 — ixi 



U7 



2 ~Vt? 



In 
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Figure 4.4: Heavy quark self-energy contributions to S for qq — > QQ. 



5,21 = 



5,12 



5,22 



= ^ln 


2tt 




a s 


In 




2tt 








71 





111 



it 



2 / ' 



— Lp — 7ii 



+- 



In 
-In 



u 



IT 



L 



Til 



(4.2.32) 



For arbitrary (3 and fixed scattering angle, we must solve for the relevant 
diagonal basis of color structure, and determine the eigenvalues. At "absolute" 
threshold, (3 = 0, we find 



pith 
1 5,11 



r^i = o, 



pth 
1 5,22 



7T 



Til 



— 1 — ixi 

2 I 2/3 



(4.2.33) 



Notice that T^ 1 is diagonalized in a basis of singlet and octet exchange in the 
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s channel, 



Csinglct Ci, 



Coctet 



'2A^ Cl+ 2' 



(4.2.34) 



with eigenvalues, 

^singlet = r^-Q , A oc tet = T^| 2 2 ■ (4.2.35) 

The general result in this s channel singlet-octet basis becomes: 



r (i.8) 

1 5,11 

r (l,8) 
1 5,21 



r (1.8) 
1 5,12 

r (l,8) 
1 5,22 



—C F (Lp + l + m), 

7T 



2a 



-^ln(^ 

7T C A 11 Ul 



7T 

4 



4 In 



Lg — 1 — 7ri 



-31n(^-ln /m2s ' 



+ L/j + iri 



(4.2.36) 



ti ) \t\U\ 

Y s is also diagonalized in this singlet-octet basis when the parton-parton cm. 
scattering angle is 9 = 90° (where Ui — ti) with eigenvalues 



Asinglet = C F (Li3 + 1 + 717 ) , 



7T 



A, 



octet 



a. 



IT 



C A 



-C F (Lp + 1 + iri) + — (L^ - In 



m 2 s 



t 2 ) + ^ 



(4.2.37) 



(4.2.38) 



It is of interest, of course, to compare the one- loop expansion of our results 
to known one-loop calculations, at the level of NLO. We may give our result 
as a function of z, since the inverse transforms are trivial. They may be found 
in terms of the Born cross section, the one-loop factoring contributions of g[ q 
and <?2 9 ^\ an d ^22^- In the DIS scheme the result is 



,(99) 



IJ 



(1) 



a. 



OBc 



7T 1 — Z 



21n(l-z) + - 
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+8 In ( j-^j - 2 - 2Re L p + 2 In ( 



(4.2.39) 



The logarithm of s//x 2 describes the evolution of the parton distributions. This 
result cannot be compared directly to the one-loop results of fl2| for arbitrary 
(3, where the singular behavior is given in terms of the variable S4, with 



S4 = (pj + k) 2 — m 2 « 2pj • . 



(4.2.40) 



where k = p a + pt, — pi — pj is the momentum carried away by the gluon. At 
partonic threshold, both S4 and (1 — z) vanish, but even for small S4, angular 
integrals over the gluon momentum with S4 held fixed are rather different than 
those with 1 — z ~ 2(pj + Pj) ■ k/s held fixed. 

Nevertheless, the cross sections become identical in the (3 —>■ limit, where 
we may make a direct comparison. Near /3 — 0, our cross section becomes 



^ J) (^,Mi,ti,s) (1) | th 

ij 



0"Bo 



71 1 — Z 



+2 In 



Am 2 



21n(l-z) + - 



(4.2.41) 



Near s = 4m 2 , we may identify 2m 2 (1 — z) — S4, and this expression becomes 
identical to the (3 —>■ limit of eq. (30) of |12|]. It is also worth noting that 
even for (3 > 0, the two cross sections remain remarkably close, differening only 
at first nonleading logarithm in the abelian (C F ) term, due to the interplay 
of angular integrals with leading singularities for the differing treatments of 
phase space. 
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As for the Drell-Yan cross section, our analysis applies not only to absolute 
threshold for the production of the heavy quarks ((3 = 0), but also to partonic 
threshold for the production of moving heavy quarks. When (5 nears unity, 
however, the anomalous dimensions themselves develop (collinear) singularities 
associated with the fragmentation of the heavy quarks, which in principle may 
be factored into nonperturbative fragmentation functions. 

Finally, we have checked that our anomalous dimension matrix for heavy 
outgoing quarks (4.2.30) reduces in the limit m — ► to the anomalous dimen- 
sion matrix for light outgoing quarks, which is || 



4.3 Applications to gg — > QQ and gg — > qq 



In this section we give the results for the anomalous dimension matrix 
when the incoming partons are gluons and the outgoing partons are heavy 
quarks 




(4.2.42) 



gij>a) + g(Pb) -> Qipi) + QiPj) ■ 



(4.3.1) 
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Figure 4.5: UV divergent one-loop contributions to S for gg — > QQ or gg. 

For the sake of completeness we also give results for the case when the outgoing 
quarks are light 

g(p a ) + g(j>b) -> q(pi) + q(Pj) • (4.3.2) 

In fig. 4.5 we show the UV divergent one-loop contributions to S for gg — > QQ 
or gg -> gg. 

Our analysis is similar to the one in the previous section. We use the 
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Figure 4.6: Eikonal rules for gg — > QQ or qq. The gluon momentum flows out 
of the eikonal lines. 

same integrals for the calculation of the The eikonal rules for incoming 
gluons are slightly modified and are given in fig. 4.6. We choose the following 
basis for the color factors: 

ci=<5<%, c 2 = d abc T^ c 3 = tf abc Tf J . (4.3.3) 

Again, the counterterms for S are the ultraviolet divergent coefficients times 
our basis color tensors: 

51 = CiZ Sj n + c 2 Z St 2i + c 3 Z s , 3i , (4.3.4) 

5 2 = c x Z s ,\i + c 2 Z s ,22 + c 3 Z s ^ 2 , (4.3.5) 

5 3 = CiZ S) i 3 + c 2 Z St 2 3 + c 3 Z S:33 . (4.3.6) 

Our results for the anomalous dimension matrix when the outgoing quarks 
are heavy are: 

r S)11 = ^jcH-^-Li-^ + lJ + ^fln^-lliiK^ + l-Tri]} 
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= o, 


rs,3i 


«s In [ Vai% 


71 \ V aj' 


^,12 


= 0, 














+L p - In 




AT 2 -A 



In ( ^) + ^ ln(v ai v bj v aj v bi ) 



2m 2 



- \n(v a v h ) + 2-ixi 



r,S,32 

rs,i3 

rs,23 
r5,33 



AN 

2 1 5,31, 

— r Si3 i 

Ts,22- 



rs,3i, 



(4.3.7) 



Again, we eliminate the gauge dependence of the heavy quarks by in- 
cluding the self-energy graphs in fig. 4.7. The contribution of the self-energy 
graphs (in the diagonal elements only) is, as before, 



C F {L % + Lj - 2). 



(4.3.8) 



In analogy to the previous section, we also have an additional —(a s /ir)CA 
in the diagonal elements. Then in terms of the invariants s, t±, and ui, the 
anomalous dimension matrix becomes 



a 



-^■[-C F (Lp + l)-C A ml 

7T 



o, 

7T 

o, 
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Figure 4.7: Heavy quark self-energy contributions to S for gg — ► QQ. 



Ts,32 
Ts,33 



M- Cp ( i(J+ l) + ^ 
iV 2 -4_ 



In 



m 2 s 



4N 

-J- 5,31, 
— ^,31, 



5,22- 



5,31 5 



+ Lg - TTt 



(4.3.9) 



At threshold the anomalous dimension matrix becomes diagonal with 
eigenvalues 



-ith 

s,n 



5,22 



5,33 



7T 

«£ 

7T 
-pth 



Til 



Tii 



C a ( 



■ 1 — Tii + 



7T2 

20, 



5,22- 



(4.3.10) 

(4.3.11) 
(4.3.12) 



We also note that the matrix is diagonalized at 9 = 90°. 

Again, it is interesting to compare the one-loop expansion of our results 



to the one-loop calculations in | l2| . In this case our calculation is complicated 
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by the fact that the color decomposition is not trivial as it was for qq. We 
have to decompose the Born cross section into three terms according to our 
color tensor basis. After some algebra and putting Cp = 4/3 and Ca = 3, our 
result becomes 



a: 



-—KggBc^isM^x) [-2881n(l - z) 



-144 In 



- 72 In 



1„2 



- 8Lp + 64 



+1281n(l-2) +64 In 



4 4 
— L - 28 - - 

9 ^ 9 



+ 28 In 



*1 Mi 

m 2 s 



where 



Sqesd(s,*i,Wi) = h — + 



1 



m 2 s 



(4.3.13) 



(4.3.14) 



U\ t\ t\U\ \ t\U\) 

The logarithms of s/fi 2 describe the evolution of the parton distributions. 
As we discussed in the previous section, our result cannot be compared 



directly to the one-loop results of [[12], but as (3 — > our expression becomes 
identical to the (3 — > limit of the sum of eqs. (36-38) in JTSJ. Even for /3 > 0, 
the two cross sections remain remarkably close. 

Finally, the anomalous dimension matrix for the case when the outgoing 
quarks are light is given by 



1 s,n 
r,s,2i 



7T I 

+C A 



In 



In ( — ) ln(z/jZ/j) + 1 — ni 



ln(v a v b ) + 1 - 7TZ 



Vab 

2 



0. 



VaiVbj 



In 

7T V^aj^&i 



5,22 



71 



+c A 



\n(v ai v bj v aj v bi ) 



2 V va 



ln(u a u b ) -ln2 + l 





AT2-4 


1^5,32 = 


4iV 


= 




1^5,23 = 




-^ F S,31, 


r5,33 = 


r^.22- 



5,31, 



(4.3.15) 

Again, we note that the matrix is diagonalized at 9 — 90°. 

Finally, we have checked that our anomalous dimension matrix for heavy 
outgoing quarks (4.3.7) reduces to the anomalous dimension matrix for light 
outgoing quarks (4.3.15) in the limit m — > 0. 



4.4 Conclusions 

We have illustrated the application of a general method for resumming 
next-to-leading logarithms at partonic threshold in QCD cross sections. We 
have given explicit results for heavy quark production through light quark 
annihilation and gluon fusion, and for light quark production through gluon 
fusion. Possible extensions include, of course, dijet and multijet production. 
We reserve estimates of the phenomenological importance of these nonleading 
terms to future work, but we hope that whether they give small contribu- 
tions or large, the method will improve the reliability of perturbative QCD 
calculations of hard scattering cross sections. 
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